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A study is reported of the same- and opposite-sign charge-dependent azimuthal correlations with 
respect to the event plane in Au+Au collisions at ys NN = 200 GeV. The charge multiplicity asym- 
metries between the up/down and left/right hemispheres relative to the event plane are utilized. The 
contributions from statistical fluctuations and detector effects were subtracted from the (co-)variance 
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of the observed charge multiplicity asymmetries. In the mid- to most-central collisions, the same- 
opposite-) sign pairs are preferentially emitted in back-to-back (aligned on the same-side) directions. 
The charge separation across the event plane, measured by the difference, A, between the like- and 
unlike-sign up/down — left/right correlations, is largest near the event plane. The difference is found 
to be proportional to the event-by-event final-state particle ellipticity (via the observed second-order 
harmonic ug bs ), where A = (1 .3± 1.4(stat)±*;°(syst)) x 1 0" 5 + (3.2 ± 0.2(stat)±£f (syst)) x lO" 3 ^ 8 
for 20-40% Au+Au collisions. The implications for the proposed chiral magnetic effect (CME) are 
discussed. 



PACS numbers: 25.75.-q, 25.75.Dw 
I. INTRODUCTION 

Relativistic heavy-ion collisions at RHIC create a hot 
and dense medium that exhibits the properties of a 
strongly coupled Quark Gluon Plasma (sQGP) [JQ. It 
is possible that chiral symmetry is restored in an sQGP. 
In addition, it has been suggested that meta-stable do- 
mains capable of undergoing topological charge changes 
can form in the sQGP, and Parity (P) and Charge con- 
jugation and Parity (CP) symmetries may be locally vio- 
lated |5l-[l2|. Kharzeev et at proposed that such a Local 
Parity Violation (LPV) can lead to the separation of pos- 
itively and negatively charged particles with respect to 
the reaction plane. This charge separation would be with 
respect to the direction defined by the axis of the classical 
magnetic field that is created by the passing ions. This 
P and CP violating process together with the magnetic 
field has been called the Chiral Magnetic Effect (CME). 

One consequence of the expected charge separation 
into pairs of back-to-back opposite-sign particles would 
be a positive two-particle azimuthal correlator, (cos(0 a + 
4>p — 2t/>rp)), of opposite-sign particle pairs and a nega- 
tive correlator of same-sign particle pairs, where 4> a and 
(j)/3 are the azimuthal angles of the two particles and t/'rp 
is the reaction plane angle [l3[ . Since the reaction plane 
angle is not known, this correlator is estimated from the 
three-particle correlator, (cos(</> Q + <j>p — 2(f> c )}, with c 
denotes the third particle, assuming that three-particle 
correlations unrelated to the reaction plane can be ne- 
glected. Previously, the STAR collaboration at the Rela- 
tivistic Heavy-Ion Collider (RHIC) measured a negative 
correlator for same-sign pairs and a small, near zero, cor- 
relator for opposite-sign pairs [2 EH- The same-sign 
result was qualitatively consistent with the expectation 
from the CME [9l-[l2|. The opposite-sign result, on the 
other hand, was inconsistent with the expectation where 
the opposite- and same-sign pair correlations should be 
equal in magnitude and opposite in sign [ij-flil]. How- 
ever, the near-zero opposite-sign result may be consis- 
tent with the CME with an additional contribution from 
in-medium interactions [l2|. More recently, the ALICE 
experiment at the Large Hadron Collider (LHC) has mea- 
sured qualitatively similar correlation signals [16|. 

It is assumed that the charge separation along the or- 
bital angular momentum axis, which is the direction of 
the magnetic field, due to the CME will induce an asym- 
metry of positively or negatively charged particle mul- 



tiplicities between the two (up and down) hemispheres 
which are separated by the reaction plane. No asym- 
metries due to the CME are expected in the left and 
right hemispheres separated by the plane normal to the 
reaction plane. In this paper, a new correlation tech- 
nique is introduced that may be sensitive to the charge 
separation that was previously investigated in STAR. 
Consistency is found between the previously published 
data [HI, [l5[ and the results from the present analysis 
when using the same charge correlator observable (see 
Appendix In this paper, a new approach is used 
to explore the charge asymmetries on an event-by-event 
basis, which extends the previous STAR measurements. 
The results obtained from d+Au and Au+Au collisions 
at v /s^ = 200 GeV at RHIC measured by the STAR 
experiment are reported [17] . 

This paper is organized as follows. The analysis 
method is described in Sec. [TTJ The data analysis tech- 
niques are described in Sec. IIIII The charge asymmetry 
correlation results are presented in Sec. IIVI The impli- 
cations of these results with respect to LPV/ CME are 
discussed in Sec. [V] Finally, the summary is presented in 
Sec. IVI1 The mathematical correspondence between the 
present charge asymmetry correlation approach and the 
previously published three- particle correlators [3 EH are 
presented in Appendix [5] Additional details about the 
data analysis are presented in Appendix [Bj The system- 
atic uncertainties of the present results are described in 
Appendix [Cj 



II. ANALYSIS METHOD 

Figure Q] schematically depicts the transverse overlap 
region of a heavy-ion collision. The event plane, denoted 
by 'EP', is reconstructed from the measured charged par- 
ticle azimuthal distributions. The event plane, 'EP', is 
not identical to the true reaction plane due to the mea- 
surement resolution. This is discussed in more detail in 
Appendix IB 61 The particle multiplicity asymmetries are 
defined on event-by-event basis via, 

^+, UD = (JV+.u- N+, U )/(N + , V + N + ,n), 

A_, UD = (iV_,u - A^ D )/(A^u + AT_ ;D ) , 

A+,lk = (V+, L - N +iR )/{N +ih + N + , R ) , and 

A_, LR = (JV_,l - N-,r)/(N-,-l + JV-,k) ■ (1) 
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Here, N +i u, jV+,d, jV+,l, and -/V+,r represent the posi- 
tively charged particle multiplicities in the up (quadrants 
I and II), down (III and IV), left (II and III), and right (I 
and IV) hemispheres as depicted in Fig. [TJ respectively. 
The same multiplicities of the negatively charged parti- 
cles are represented by N- \j, N^^, 2V l, and N-^ R . 




FIG. 1: Schematic depiction of the transverse overlap region 
of a heavy-ion collision. The event plane (EP) direction is 
reconstructed from final state particle momentum space. Four 
quadrants are defined and are labeled: Up = I+II, Down = 
III+IV, Left = 11+111, Right = I+IV. 

The topological charge-signs are expected to be ran- 
dom from one metastable domain to another in a single 
event and in different events [||-[l2j]. Also, the reaction 
plane obtained from the second harmonic of the particle 
azimuthal distributions cannot distinguish up from down. 
Thus, while the magnitude of the up-down (UD) multi- 
plicity asymmetry becomes larger due to the LPV/CME, 
its sign is random. As a result, the average asymme- 
tries remain zero, but the distributions of j4±,ud, where 
^4±,ud is used to denote A +! ud and A^ud, would be 
wider than those of A±^ R , (i.e. A + .lr and A_.lr)- In 
other words, the variances (A± UD ) should be larger than 
the variances (Aj_ LR ). Therefore, the variances of the 
charge multiplicity asymmetries is of interest here. 

The covariance of the charge multiplicity asymme- 
tries, (A + A-) VD = (A+^udA- UD ) and (A+A-) LR = 
(A-^lr^-.lr), are also studied. The covariance mea- 
sure is analogous to the traditional parity violation mea- 
sures as follows. The multiplicity asymmetry of one 
charge sign, e.g. A +! ud, either positive or negative, can 
be viewed as to define preferentially the "parity-axis" 
direction, the combined effect of the orbital angular mo- 
mentum direction and the topological charge sign. The 
covariance, (A + A_)ud, is then a measure of A^ t \ju with 
respect to this "parity-axis." 

The charge asymmetry correlations are, however, 
parity-even and subject to physics backgrounds similar 



to those in the charge correlator measurement. These 
physics backgrounds can be assessed by the left-right 
(LR) asymmetry correlations, (A^) and (A + A-)i, R , to 
which the LPV/CME does not contribute. The LR 
charge asymmetry correlations could thus serve as the 
null-result reference. However, as will be discussed in 
Sec. IV B[ the physics backgrounds to the UD and LR 
measurements may be different. 

In the present analysis, the charge multiplicity asym- 
metries, ^4±,ud and j4±,lr> are computed event-by-event. 
The variances of these quantities, (A 2 ^ y D ) and {A^_ ^ R ), 
are obtained and the statistical fluctuations and detector 
effects are subtracted from these, via 



S(Al LR ) 



— (^Iud) — (^±,UD,stat+dct) 
= (^±,Lr) — (^±,LR,stat+det) 



(2) 



The details on the contributions from the statistical and 
detector effects are presented in Appendix IB 41 A com- 
parison of 6(A+) and £>(A?_) is described in Appendix lB 51 
and the two quantities are consistent. Thus, average dy- 
namical variances, 

5{A 2 VD ) = (<5« UD ) + <5(A 2 _, UD ))/2, 
5{Al K ) = (S(A 2 +tLR ) + S(A^ LR ))/2, (3) 

are presented. The covariances, (^4+j4_)ud and 
(^4 + A_)lr, are also presented. The statistical fluctua- 
tions and detector effects on the covariances are consis- 
tent with zero. Nevertheless, they are removed from the 
covariances: 

5(A+A-) VD = (A +J 4_) UD - (A + ^_) UDiStat+dot , 
S(A + A-) LR = (A + A-) LR - (A+A-) LR , stat+dct . (4) 

The differences between the UD and LR measurements 
which may be directly sensitive to the CME will be re- 
ported. Namely, 



A(A 2 ) = S(A 2 VB ) - S(A 2 LR ) , 
A(A + A_) = 5(A + A_) VU - 8{A + A_) 



LR • 



(5) 



are the same as (A 2 ] - d ) — (A 2 jR i 



The A(A 2 ) and A(A + A- 

and ( j 4 + A_)ud — (A+A-)i, R , resepectively, because the 
statistical fluctuation and detector effects cancel in the 
differences. The difference UD— LR correlations between 
same- and opposite-sign charges, 



A = A(A 2 ) - A(A+A_> 



(6) 



will also be reported, which may quantify the charge sep- 
aration effects. 

In addition to the asymmetries between hemispheres, 
the asymmetries between azimuthal wedges of smaller 
sizes will be studied as depicted in Figure [U(a). The out- 
of-plane asymmetries between a wedge of size 2A0 W at 
4> w = 90° and the same-size wedge at <f> w = 270° will 
be explored. This involves the counting of the particle 
multiplicity within an azimuth range relative to the EP 
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between 90° — A0 W and 90° + A</> w , and the multiplicity 
within 270°-A</> w and 270°+A</> w . Similarly for in-plane 
asymmetries, the particle multiplicities within 0° — A</> w 
and 0° + Acj) w , and within 180° - A0 W and 180° + A0 W 
are counted. 




FIG. 2: Schematic view of the transverse overlap region and 
the back-to-back wedges (azimuthal regions) where the charge 
asymmetries are computed. The event plane (EP) direction is 
reconstructed from final state particle momentum space, (a) 
Configuration for the study of the wedge size dependence, (b) 
Configuration for the wedge location dependence. 

The charge asymmetry correlations will also be cal- 
culated within a back-to-back pair of wedges at spe- 
cific azimuthal locations with respect to the EP. Fig- 
ure [SJb) shows the schematic configuration of the back- 
to-back wedges with size 2A0 W at location <j) w . For 
these studies, the notations used for the asymmetry cor- 
relation variables are refined as follows. The variable 
^±,0w±A<ji w stands for the positive/negative particle mul- 
tiplicity asymmetry between the wedge (f> w ± A</> w and its 
opposite-side partner wedge (ir + <f) w ) ± A</> w . In this no- 
tation, the variables A-^ud and ^4±,lr are equivalent to 
-4±,90°±90° and A ± . o± 9 o°, respectively. 

A disadvantage of using smaller, non-hemispheric, 
ranges is the fact that the statistical fluctuation and 
detector effects no longer cancel between variances 
from out-of-plane wedges and those from in-plane 
wedges. The statistical fluctuation contributions to 

( A ±,0 w ±A0 w ,stat+det) and (^+^-)^±A^,,stat+det are ob- 
tained in the same way as described in Appendix IB 41 
These contributions are subtracted to obtain the dynam- 
ical correlation: 



Similarly, the average 



(^±,0 W ±A0 W ) — 

,0 w ±A0 w ,stat+det/ ' 
(A+A-) ( p w ±A<p w - 

w ±A0 w ,stat+det ■ 



(7) 



,±A0„ 



-<*«0 w ±A0 w ))/2, 



(8) 



and the differences between the UD and LR correlations 
will be reported. 

The charge multiplicity asymmetry correlation quan- 
tities in the present study can be related to the three- 
particle charge-dependent azimuthal correlators [TJ], UH 
as shown in Appendix DO The present analysis differs 
from that in Refs. [HI, Tla ] in that the reported corre- 
lations represent the contributions from the entire cor- 
relation structure (all multipoles), whereas those in the 
previous studies focused on the lowest-order azimuthal 
multipole only. 



III. DATA ANALYSIS 

The data used in this analysis were taken by the STAR 
experiment [l8j at RHIC at the nucleon-nucleon center 
of mass energy of ^/s NN = 200 GeV. The minimum-bias 
and central triggered Au+Au data were from RHIC Run 
IV in the year 2004. The reference minimum-bias d+Au 
data used for comparison were from RHIC Run III in 
2003. The Au+Au results will also be compared to data 
taken from Run VII (2007). The Run VII data were used 
for the study of the asymmetry correlations with respect 
to the first harmonic event plane. 

The details of the STAR experiment can be found in 
Ref. (l8| . The minimum-bias triggers for Au+Au and 
d+Au collisions were provided by the Central Trigger 
Barrel [l9j] and the Zero Degree Calorimeters (ZDC) [20j | . 
A total of 8.8 million Au+Au events with centrality rang- 
ing from 0-80% (Run IV) and a total of 8.9 million d+Au 
events (Run III) were used in this analysis. The Run VII 
data used for comparison and for the ZDC first harmonic 
event plane study consists of 70 million minimum-bias 
Au+Au events. 

The Au+Au collision centrality is defined according to 
the measured charged particle multiplicity in the main 
Time Projection Chamber (TPC) within the pseudo- 
rapidity range |ry| < 0.5 [2l|. Results will be presented as 
a function of centrality in terms of the number of partic- 
ipant nucleons, -/V par t, which was obtained from a Monte 
Carlo Glauber calculation (22|. 

The main detector used for this analysis was the Time 
Projection Chamber (TPC) [H H. The TPC is sur- 
rounded by a solenoidal magnet providing a uniform mag- 
netic field of 0.5 tesla along the beam direction. Particle 
tracks were reconstructed by the TPC. The primary col- 
lision vertex was reconstructed using TPC tracks passing 
various quality cuts. Events with a primary vertex loca- 
tion within ±30 cm of the geometric center of the TPC 
along the beam axis were used in the analysis. 

In the present asymmetry calculations and TPC event- 
plane construction, only those "primary" tracks extrap- 
olated to within 2 cm of the primary vertex were used. 
The tracks were required to have at least 20 (out of a 
maximum of 45) hits in the TPC used in track recon- 
struction. The ratio of the number of hits used in the 
track reconstruction to the maximum possible number of 
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hits for a given track was required to be greater than 0.51 
to eliminate multiple track segments reconstructed from 
a single particle trajectory. These cuts were varied to 
assess the systematic uncertainties in the present results 
which are discussed in Appendix [C] 

The second Fourier harmonic in azimuth was used 
to determine the event-plane angle ipEP 25| from the 
TPC-reconstructed tracks. The event plane was recon- 
structed for both Au+Au and d+Au collisions, and does 
not necessarily correspond to a specific plane in configu- 
ration geometry. The transverse momentum, px, range 
of the particles used to determine the event plane was 
0.15 < pr < 2 GeV/c. The low-px cut-off was imposed 
by the magnetic field strength and the TPC inner ra- 
dius. The pT-weight method |25| was used for the event 
plane reconstruction as it gives a better event plane res- 
olution than no p^-weight in the presence of the stronger 
anisotropy at larger values of px ■ The event plane recon- 
struction was done in two different rj ranges: —1 < i] < 
and < rj < 1 (additional details below) . The slight non- 
uniformities of the efficiency and acceptance in azimuth 
due to the TPC sector boundaries was corrected for in 
the event-plane construction by using ^-dependent effi- 
ciencies (see Appendix IB II) . The azimuthal angle of the 
EP constructed by the second harmonic ranges from to 
7r. In half of the events chosen randomly, tt was added 
to the reconstructed EP azimuthal angle so that the re- 
sulting EP azimuthal angle ranges from to 2ir. The 
first-order harmonic event plane was also measured in- 
dependently using the ZDC Shower Maximum Detector 
(SMD). The ZDC-SMD event-plane analysis is described 
in Ref. [H H, El • 

The particle azimuthal angle <j> relative to EP was 
properly folded with -0ep so that <f> — i/'ep was also in 
the range between and 2tt. A particle is assigned 
to the 'up' hemisphere if < <f> — t/>ep < 7r, 'down' 
hemisphere if tt < <f> — ?/>ep < 2-7T, 'left' hemisphere 
if 7r/2 < <j> — V-'ep < 37r/2, and 'right' hemisphere if 
3tt/2 < - Vep < 2tt or < - Vep < n/2. The 
number of particles, weighted by the efficiency correction 
factor (described in Appendix IB 3p . in the upper, lower, 
left, and right hemispheres was counted. The asymme- 
tries from the corrected particle multiplicities were calcu- 
lated. Separate calculations of the asymmetries (i) using 
particles within < rj < 1 with the EP constructed from 
— 1 < i) < 0, and (ii) using particles within —1 < r\ < 
with EP constructed from < i] < 1, were performed. 
These two results were consistent (see Appendix IB 5[) . so 
their average is discussed. 



IV. RESULTS 



random from event to event. The data indeed show zero 
single asymmetries within the statistical errors. 

Figure [3] shows the dynamical variances, S(A 2 J - D ) (solid 
squares) and S(A 2 jR ) (hollow squares), and covariances, 
5(A + A-)uy) (solid circles) and 5( j 4 + ^4_)lr (hollow cir- 
cles), as a function of -/Vp art . Since two-particle correla- 
tion measures are typically diluted by a multiplicity fac- 
tor, the dynamical fluctuation quantities are multiplied 
by A^part to reveal the magnitudes over the entire cen- 
trality range. The UD and LR quantities are different 
within all centralities ranges except the most peripheral 
and most central collisions. If the EP were random and 
unrelated to the reaction plane, then the UD and LR ob- 
servables would be the same within the statistical uncer- 
tainties. This is not the clearly shown in Fig. [3] 
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FIG. 3: (Color online) Centrality dependences of the charge 
asymmetry dynamical correlations, S{A 2 ), and the positive 
and negative charge asymmetry correlations, 5{A+A-). The 
asymmetries are calculated between hemispheres separated 
by the event plane (UD) and between those separated by the 
plane perpendicular to the event plane (LR). The asymme- 
try correlations are multiplied by the number of participants 
A'part- The error bars are statistical only. The upper (lower) 
shaded band shows half of the systematic uncertainty in the 
5(A+A-} (8{A 2 )); the larger of the UD and LR systematic 
uncertainties is drawn. The stars and triangles depict the 
d+Au results. 



A. Charge Asymmetry Correlations 

The single asymmetry quantities (A+ ud); (A— ^ue)), 
(j4 +; lr), (-<4_.lr) are, by definition, zero because the po- 
sitions of the up (left) and down (right) hemispheres are 



A positive 5 (A 2 ) indicates a broadening of the asym- 
metry distributions of A + and A_ due to dynamical pro- 
cesses, whereas a negative 8 (A 2 ) indicates narrowing of 
the distributions. In peripheral collisions, both S(A 2 J - D ) 
and <5(^4lr) are positive, suggesting that the same-sign 
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particles within one unit of pseudo-rapidity \r)\ < 1 are 
more likely emitted in the same direction. This small- 
angle correlation is stronger in the up-down hemispheres 
than in the left-right hemispheres. The small-angle cor- 
relation becomes weaker when the collisions are more 
central. In fact, in more central collisions, the 5(A 2 ] - D ) 
and ^(^lr) become negative, i.e., the same-sign charge 
pairs are preferentially emitted back-to-back in those col- 
lisions. 

The correlations between A + and A-, both UD and 
LR, are large and positive implying strong correlations. 
The correlation is on the order of ~ 10 -3 , suggesting that 
the correlated asymmetry is as large as a few percent. 

Figure |3] also shows the asymmetry correlations in 
d+Au collisions. The d+Au data lie at the endpoint of 
the Au+Au curve and are consistent with an extrapola- 
tion of that trend. 

The data in Fig. [3] seems to indicate the following pic- 
ture. In d+Au and peripheral Au+Au collisions, the par- 
ticles within one unit of pseudo-rapidity are preferentially 
emitted in the same direction, whether they are the same 
or opposite charge signs. The magnitude of the small- 
angle correlation is, however, stronger in the opposite- 
than in the same-sign pairs, and is stronger in the out-of- 
plane than in the in-plane direction. In medium-central 
to central collisions, while the opposite-sign pairs are still 
preferentially aligned in the same direction and more so 
than in peripheral collisions, the same-sign pairs are pref- 
erentially back-to-back. The small-angle correlation be- 
tween the opposite-sign pairs is always stronger out-of- 
plane than in-plane. The tendency of back-to-back emis- 
sion of same-sign particles is weaker in the out-of-plane 
than in the in-plane direction. 

In order to investigate the possible contributions from 
the LPV/CME, the difference between the UD and LR 
asymmetry correlations, A(A 2 ) and A(A + A_), was stud- 
ied. The contributions from the detector effects and sys- 
tematics are largely canceled in these differences (see Ap- 
pendix [Cj. Figure H shows A(A 2 ) and A(A + A_) as a 
function of the centrality. The upper shaded area shows 
the systematic uncertainty in A(A + A-) and the lower 
shaded area shows that on A(A 2 ). Also shown as the 
lines are the A(A 2 ) and A(A + A-) values which would 
be expected with a perfect event-plane resolution, which 
is calculated assuming the linear extrapolations shown in 
Fig. [24] (right panel). As noted in Appendix \A~[ this lin- 
ear extrapolation would be correct if high-order harmonic 
terms in Eq. (|A2I) are negligible. 

Figure 2] shows that A(A 2 ) is larger than zero, i.e. 
5(A 2 j D ) > S(Af iR ), in collisions at all centralities. The 
UD asymmetry distribution is always broader than the 
LR one. This implies that there are more small-angle 
same-sign pairs in the out-of-plane direction than in the 
in-plane direction. Equivalently, there are more back- 
to-back pairs in-plane than out-of-plane. The UD— LR 
difference in 5(A + A-) is small relative to the correla- 
tions themselves. This indicates that the majority of 
the strong correlations between opposite-sign pairs is 
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FIG. 4: (Color online) The correlation differences A(A 2 ) = 
6{Al D }-5{Al R ) and A(A+A-} = 5(A + A-} VD -S(A+A-} LR , 
scaled by the number of participants iVp a rt, as a function of 
A'part- The error bars are statistical, and the systematic un- 
certainties are shown in the shaded bands (upper band for 
A{A+A-} and lower band for A{A 2 }). Also shown as the lines 
are the linear-extrapolated values of A{A 2 } and A{A+A-} 
corresponding to a perfect event-plane resolution. The star 
and triangle depict the d+Au results. 



unrelated to the reaction plane. On the other hand, 
A(A + A-) is also larger than zero, i.e. <5(A + A_)ud > 
6(A + A-)L-[i in all centralities except peripheral colli- 
sions. The opposite-sign particle pairs are more strongly 
emitted in the same direction out-of-plane compared to 
in-plane. 

The general trend of A(A 2 ) with centrality is as fol- 
lows. It increases with centrality, reaching a maximum 
in mid-central collisions, and then decreases with increas- 
ing centrality. The trend for A(A + AJ) is similar, except 
that it is slighlty negative in peripheral collisions. The 
negative values of A(^4 + ^4_) in peripheral collisions (and 
also in d+Au collisions) are likely due to "non-flow" (e.g. 
di-jets). The decrease for more central collisions may be 
due to the experimental event-plane resolution because 
the difference should disappear in zero impact parameter 
collisions where the reaction plane is undefined. 

Figure O shows the px dependence of the charge asym- 
metry correlations in 20-40% and 0-20% central Au+Au 
collisions. The values of 8 (A 2 ) are positive at low px, and 
decrease sharply with increasing px up to about 1 GeV/c. 
In this centrality bin, the 8(A 2 jB ) values become negative 
and S(A 2 j- D ) is approximately zero for px > 1 GeV/c. 
This indicates that low-p^ pairs are emitted preferen- 
tially in the same direction. The back-to-back emission 
of same-sign pairs increases with increasing px- On the 
other hand, 5(A + AJ) remains relatively constant at low 
Px up to roughly 1 GeV/c, and then increases sharply 
with increasing px- The aligned emission of opposite- 
sign pairs in the same direction increases strongly with 
px above 1 GeV /c. The observed features in px are qual- 
itatively similar for other centrality bins. It is worthwhile 
to note that the charge asymmetry correlations in each 
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Pt bin, shown in Fig. [5l are calculated solely from the 
particles within that pr bin. On the other hand, the 
charge asymmetry correlations as a function of central- 
ity, shown in Fig. [31 are calculated from all particle pairs 
with pr < 2 GeV/c. As such, the results shown in Fig. [3] 
for a given centrality bin cannot be trivially obtained 
from those in Fig. [5] for the same centrality bin. 
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FIG. 5: (Color online) The pr dependence of the charge asym- 
metry dynamical correlations, 5 {A 2 }, and the positive and 
negative charge asymmetry correlations, 5(A+A-). The data 
are from 20-40% central Au+Au collisions. The asymmetries 
are calculated between hemispheres separated by the event 
plane (UD) and between those separated by the plane perpen- 
dicular to the event plane (LR). The error bars are statistical 
only. 

Figure [6] shows the UD— LR correlations as a function 
of pt for 20-40% central Au+Au collisions. The val- 
ues of A(yl 2 ) and A(A + A-) increase with pr- There 
is qualitatively no difference in the pt dependence be- 
tween the same- and opposite-sign correlations. On the 
other hand, it is generally expected that the UD— LR dif- 
ference should be most prominent at low pt if the CME 
is responsible. Such a low-pr feature is not observed in 
these data. This was als o q ualitatively observed in the 
three-particle correlators [lj, [l5| ■ 



B. Dependence on Event- by- Event Anisotropies 

There could be many physics mechanisms con- 
tributing to the event-plane dependent charge correla- 
tions. Voloshin estimated that resonance decays with 
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FIG. 6: (Color online) The correlation differences A{A 2 ) = 
^ud)-<5<^lr) and A<A+A_> = 6(A + A-} VD -6(A + A-} LR 
as a function of pr- The data are from 20-40% central (upper) 
and 0-20% central (lower) Au+Au collisions. The error bars 
are statistical only. 



anisotropies were insignificant to charge dependent cor- 
relations relative to the reaction plane Wang sug- 
gested that general cluster particle correlations with 
anisotropies could generate a sizable difference between 
the in-plane and out-of-plane particle correlations [28| . 
Pratt et al. argued that momentum conservation and lo- 
cal charge conservation together with elliptic flow could 
yield event-plane dependent correlations that differ be- 
tween same- and opposite-sign pairs J29(. A path- length 
dependent jet-quenching effect (30M35j j could also con- 
tribute [33. To test these ideas experimentally, the 
values of 8(A 2 ), 5(A + A-}, A(A 2 ), and A(A +J 4_) were 
studied as a function of the azimuthal elliptic anisotropy 
of high-pT and low-pT particles. The high-p^ anisotropy 
may be most sensitive to the jet-quenching effect, while 
the low-pT anisotropy characterizes the bulk event shape. 
The event elliptic anisotropy is computed via (cos2(</> — 
V'ep)) at low pt and high pt- For low pr, only the par- 
ticles within one half of the TPC (also used in the asym- 
metry measurements) were used to compute the elliptic 
anisotropy, while the angle -0ep was reconstructed using 
the particles in the other half of the TPC. The variable 
i>£ bs is used to stand for the \ow~pt anisotropy and is de- 
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FIG. 7: (Color online) The charge asymmetry correlations 5(A 2 } (left panels) and S(A+A-) (center panels), and correlation 
differences A{A 2 } — 8{A 2 J - D ) —S{A 2 j - R ) and A{A+A-} = 5{A+A-)\j-d — <5(A+ j 4_)lr (right panels), as a function of the azimuthal 
elliptic anisotropy of high-pT (pr > 2 GeV/c) particles (upper panels) and low-pr (pr < 2 GeV/c) particles (lower panels). 
Data are from 20-40% Au+Au collisions. The particle pr range of 0.15 < pr < 2 GeV /c is used for both EP construction and 
asymmetry calculation. Error bars are statistical. 



fined via v% hs = (cos 2((f>(p T ) - ^ep))|o.15< P t<2 GcV/c- 
For high pt, particles with px > 2 GeV/c from the 
entire TPC (\r]\ < 1) were used to increase the statis- 
tics. The variable ^2pT>2GeV/c ^ s usec ^ *° s t an d for the 
high-px anisotropy and is defined via ^2p T >2GeV/c = 

(cOs2((/)(p T ) - '0Ep))|p T >2 GoV/c- 

Figure [Jj (upper panels) shows the asymmetry cor- 
relation results in 20-40% central Au+Au collisions as 
a function of v^ T>2 GeV / c - Over the large range in 
t, 2pT>2GcV/ci relatively small variations are observed. 
This may indicate that path-length dependent jet- 
quenching does not have a significant impact on the 
charge asymmetry correlations. This is consistent with 
the theoretical study described in Ref. |36| . 

Figure [JJ (lower panels) shows the asymmetry correla- 
tion results as a function of v% hs . Significant changes are 
observed in the variances. The values ofd(Al R ) decrease 
with increasing v% hs while the values of 5(Ay D ) increase. 
This results in a strong increase in the difference A (A 2 ) 
with increasing i>2 bs . The variations in the covariances 
are significantly weaker. This results in a weaker depen- 
dence of A(A + A-) on v% hs . However, the change appears 
to be in the opposite direction, decreasing with increasing 

obs 

v 2 ■ 

It is interesting to note that the A (A 2 ) and A(A+A-) 
cross at v% hs and the crossing point is at a positive 
value. For v% hs = 0, no difference is apparent between the 



same- and opposite-sign pair correlations. Of course, the 
average v% hs in these data is nonzero, but positive. As 
a result, the asymmetry variance of all events is larger 
than the covariance. As suggested in Appendix IB 61 it 
is possible that, in events with negative V2 hs , the recon- 
structed EP does not reflect the true reaction plane, per- 
haps being orthogonal to rather than aligned with the 
reaction plane. This would mean that the UD and LR 
are flipped for those events with a significantly negative 
^obs rp ne magnitudes of the asymmetry correlations de- 
pend on centrality, but their qualitative features versus 
u 2p T >2GoV/c an d v 2 hs are similar for different centrali- 
ties. 

The measurements shown in the lower panels of Fig. [7J 
were also performed with a random event plane. The 
results with the random event plane are similar to those 
shown in Fig. [Jj They have the same dependence on 
V2 hs , but the same- and opposite-charge results in the 
right lower panel cross at zero intercept. 



C. Dependence on Wedge Size 

Above, the charge multiplicity asymmetry correlations 
between hemispheres have been described. One advan- 
tage of these asymmetries is that they are calculated 
from the same set of particles, and are only divided ei- 
ther UD or LR. The statistical fluctuations and detector 
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effects cancel in the difference between UD and LR, so 
A(A 2 ) ee 5{A* D ) - d{Atx) = (Alv) - (At R ). However, 
measurements of multiplicity fluctuations within hemi- 
spheres are not sensitive to possibly smaller scale angular 
structures of the charge separation. For example, the cor- 
related charged particle pairs from the CME that were 
initially aligned with the total angular momentum di- 
rection may or may not remain aligned (or preferentially 
aligned) in the same direction [l2| . In order to investigate 
the angular structure of the charge separation, the charge 
multiplicity asymmetry measurements were restricted to 
azimuthal ranges ("wedges"), A</> w , that are smaller than 
7r (hemispheres), allowing the study of the charge sepa- 
ration as a function of the wedge size. 

Figure [3] (upper panel) shows the asymmetry correla- 
tions, S(Al 0O±A(}> J, S(Al 0±A<j) J, (5(A + A_) 9O o ± A0 w , and 
(5(^4+A_)o°±A0 w , versus the wedge azimuthal size (see 
Fig. [2fa)). Both the covariances, S(A + A^) 90 °±A4, V , and 
<5(yl + A_)o°±A0 w , increase with decreasing wedge size, 
A0 W . This suggests that the major contribution to the 
opposite-sign charge correlation is local. The variances, 
(5(^4g O o ± A0 ) an d <5(^o°±A0 )) increase with decreasing 
A0w- 



asymmetry correlations, A(A A(j) ^) = <5(^4g O °±A0 w ) — 
$(Alo± A( pJ and A(A + /1_)a0 w = 5{A + A_)9 o± A<t>w - 
5(A + yl_)o°±A0 w , between the out-of-plane and in-plane 
directions. Both A{A A(tl J) and A(A + A^)a4>„ increase 
with decreasing wedge size and have qualitatively similar 
trends. The difference between the two seems to diminish 
with decreasing wedge size. 

In the above, the focus has been on the difference 
between the in-plane and out-of-plane wedges. In the 
following, the charge asymmetry correlations in fixed- 
size back-to-back wedges, as a function of the wedge 
azimuth relative to the event plane, are discussed (c/. 
Fig- Hlb)). Figure IH1 shows the asymmetry correlations 
between 30°-wide back-to-back wedges, §(A% ± 150 ) and 
5(A + j4_)0 w ±i5o , versus the wedge azimuthal location 
4> w relative to the event plane. The data are from 20- 
40% central Au+Au collisions. The asymmetry cor- 
relations increase from in-plane to out-of-plane, as ex- 
pected, for both same-sign and opposite-sign charges. 
The EP-independent part of the correlations is stronger 
in opposite-sign charges. The dependence appears to fol- 
low the characteristic behavior of cos(20 w ). 
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FIG. 8: (Color online) The wedge size dependences of charge 
multiplicity asymmetry correlations (upper panel), and their 
differences between out-of-plane and in-plane, A(A A w ) and 
A(A + 4_) A W (lower panel) for 20-40% central Au+Au colli- 
sions. The error bars are statistical only. 

Figure [5] (lower panel) shows the difference in the 



FIG. 9: (Color online) Charge multiplicity asymmetry cor- 
relations as a function of the wedge location, <^) w , in 20-40% 
central Au+Au collisions. The wedge size is 30°. The curves 
are the characteristic cos(2<j!> w ) to guide the eye. The error 
bars are statistical only. 



V. DISCUSSIONS 

These measurements were motivated by the 
LPV/CME. The LPV/CME produces quark charge 
separations along the system's magnetic field axis |a-[lO| . 
This may result in a charge separation reflected by the 
final state hadrons. A negative correlator for same-sign 
pairs was observed [13, [H| . This was qualitatively con- 
sistent with the expectation from the local strong parity 
violation However, a small, close to zero, corre- 

lator for opposite-sign pairs was also observed 0, EU ■ 
The result appeared to be inconsistent with the naive 
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CME expectation above. To explain the preliminary 
version of the results in Ref . [H], [H[ , Kharzeev et ah [Hfl 
suggested that the initial correlation among the quarks 
from CME would not survive the subsequent dynamic 
evolution to the final state hadrons, except those emitted 
from the surface of the collision zone. The back-to-back 
correlations of opposite-sign pairs from CME would 
be lost because at least one quark of a pair would be 
affected by interactions with the evolving medium. 
Such a medium interaction scenario could qualitati vely 
explain the observed opposite-sign correlator [l2l H3. Il5j. 
More recent measurements by STAR [33| may indicate 
that a parity conserving background is more likely to 
explain the suppression rather than the medium induced 
back-to-back suppression previously supposed [HI, [H| • 

The LPV/CME should produce more same-sign charge 
pairs in the up-down hemispheres, yielding wider (larger 
variance) asymmetry distributions of both positive and 
negative particle multiplicities. Therefore, the up-down 
dynamical variance S(A 2 j- D ) should be larger than the 
left-right dynamical variance S(Al R ). The LPV/CME 
should also produce an anticorrelation between multiplic- 
ity asymmetries of positive particles and negative parti- 
cles in the up-down rather than in the left-right hemi- 
sphere. This should result in <5(A +j 4_)ud being smaller 
than 6(A + A-) L r. 

The present results show 5(A 2 JD ) > S{Af iR ) in Au-f Au 
collisions at all centralities at 200 GeV (see e.g. Fig. 0}. 
The UD asymmetry distribution is wider than the LR 
asymmetry distribution. This is qualitatively consistent 
with the CME, although the same-sign particles are pref- 
erentially back-to-back in medium-central to central col- 
lisions regardless of whether they are oriented in-plane 
or out-of-plane. However, (5(A +j 4_)ud > 8(A+A-)i,-r in 
all centralities except peripheral collisions. The opposite- 
sign particle pairs are more strongly emitted in the same 
hemisphere in the UD than LR direction. This seems con- 
trary to the expectations from the CME alone. In d+Au 
and very peripheral Au+Au collisions, S(A + A^)u^, < 
<5(A + j4_)lr- In these small multiplicity collisions, there 
might exist autocorrelations between particles in the 
t] < and 1] > regions. One such constraint may be 
due to "non-flow," e.g. di-jets, causing the reconstructed 
event plane to align preferentially with the plane con- 
taining the two jets. In this case, the left and right hemi- 
spheres are the two jet-hemispheres, which have large 
fluctuations in multiplicity and result in large LR asym- 
metries. The up and down hemispheres, on the other 
hand, are a roughly symmetric division of the event with 
respect to the di-jet axis, resulting in smaller asymme- 
tries. 

There is little theoretical guidance regarding the quan- 
titative magnitude of the charge asymmetry expected 
from the CME. An order of magnitude estimate suggests 
that the c harg e asymmetry could be on the order of a few 
percent [9l4llj , charge asymmetry correlations on the or- 
der 10~ 4 -10~ 3 . The medium attenuation effects likely 
reduce the asymmetry correlations by an order of mag- 



nitude [12j. Recently, it has been argued that charge 
asymmetries can arise from strong magnetic fields and 
well-known QCD processes without invoking LPV [38| . 
However, the estimated magnitude of the charge asym- 
metry correlations from those processes is orders of mag- 
nitude smaller [38j . 

It is worthwhile to point out that the present charge 
multiplicity correlation observables are connected to the 
three-particle correlators in Refs. [3 EH, as described 
in detail in Appendix [A] The correlation observables re- 
ported in this paper include the entire correlation struc- 
ture while the three-particle correlators focussed on the 
lowest-order azimuth multipole only. The comparison of 
the two observables suggests that the higher order multi- 
poles may be significant in the opposite-sign charge cor- 
relations, but they are insignificant in the same-sign cor- 
relations. When analyzed using the same three-particle 
correlator observable, the present data are consistent 
with those in Refs. [II El. 

It is also worthwhile to point out that the "modulated 
sign correlations" recently reported in Ref. [13] are more 
closely related to the present charge multiplicity corre- 
lation observables than the three-particle correlators in 
Refs. [l3. Hlf. In fact, the opposite-sign modulated sign 
correlation in Ref. (37| is identical to the present opposite- 
sign multiplicity correlation observable A(A+A-) except 
for a constant multiplicative factor. The opposite-sign 
modulated sign correlation from the Run VII data in 
Ref. [33] is consistent within errors with the Run IV 
A(A + A-) data reported here. 



A. Is the Charge Separation In- or Out-of-Plane? 

Although both the same- and opposite-sign UD— LR 
correlation results are positive in non-peripheral colli- 
sions (see Fig. 2]), the same-sign result, A(A 2 ), is larger 
than the opposite-sign result, A(A+A-). It is possi- 
ble that the various backgrounds may produce correla- 
tions that fall in between - being equal for the different 
sign pairs (i.e. the physics backgrounds give zero charge 
separation). Then, the different results for same- and 
opposite-sign charges, now having different signs once 
subtracting the common background, would be consis- 
tent with CME. No medium effects such as those de- 
scribed in Ref. [l2| would be needed. One may simply 
refer to the larger UD— LR of same-sign than opposite- 
sign correlation, A = A(A 2 } — A(A + A-) > 0, as "charge 
separation" across the event plane. 

If the CME results in same-sign pairs in the final 
state that are still preferentially directed along the mag- 
netic field axis [l2T| . or are of the form a\ sin(0 Q + (ftp — 
2'0rp) [l34T5j . then the difference should depend on the 
wedge size. The smaller the wedge size, the larger the 
CME effect on the difference A(A0 w ) = A(A\ lf) J - 
A(A + ^4_)a0 w - However, it is also possible that, due to fi- 
nal state interactions, the charge separation effect across 
the reaction plane is no longer preferentially directed in 
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the orbital angular momentum direction, so that A(A(f> w ) 
may not increase with decreasing A(f> w . Nevertheless, it 
is interesting to examine A(A0 w ) as a function of A0 W 
which may reveal the angular distribution of the charge 
separation effect. 

Figure [TU] shows the values of A(A0 w ) versus A0 W . 
The charge separation effect, A(A0 w ), decreases with 
decreasing wedge size, A</> w . It appears that the differ- 
ence between same- and opposite-sign correlation is di- 
minished for small wedge sizes. The largest difference be- 
tween the same- and opposite-sign pairs is obtained when 
whole hemispheres are used. This suggests that the ef- 
fect of charge separation across the event plane happens 
in the vicinity of the in-plane direction rather than out- 
of-plane. 



x10 




20 40 60 80 
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FIG. 10: The wedge size dependence of the difference between 
the same-sign and opposite-sign, A(yl^ w ) — A(A+A-)a4,„, 
shown in Fig. [8] The error bars are statistical only. 



The cartoon in Fig. [TT] depicts a map of the signal 
charges consistent with the present results. The sup- 
posed common background used in this discussion, which 
falls in-between the same- and opposite-sign measure- 
ments, is excluded from the cartoon. There are pref- 
erentially more back-to-back same-sign pairs along the 
in-plane direction. The positive pairs preferentially oc- 
cupy one hemisphere either above or below the reaction 
plane and the negative pairs preferentially occupy the 
opposite hemisphere. If those hadron pairs are the re- 
sult of the CME from the initial quarks, then the data 
would suggest that those quarks initially moving perpen- 
dicular to the reaction plane (along the magnetic field 
direction) [l2j have been deflected (and hadronize) to- 
ward the reaction plane direction. The same-sign cor- 
relator measure that was previously reported does not 
distinguish between deflected pairs and the initial pairs 
without any deflection [3!| . The opposite-sign correlator 
would be close to zero for the configuration depicted in 
the cartoon. 




- Px(EP) 



FIG. 11: Schematic depiction of the charge pairs responsi- 
ble for the observed wedge-size dependence of the difference 
between the same- and opposite-sign UD— LR correlations. 
Note that the same- and opposite-sign pairs that would yield 
the assumed common physics background falling between the 
same- and opposite-sign UD— LR measurements are excluded 
from this cartoon. 



Is the Charge Separation a CME Signal or 
Background? 



The charge separation effects described above could be 
due to the CME if the physics backgrounds for the same- 
and opposite-sign pair correlations are similar and fall in 
between the two, although the charge separation seems 
to happen in the vicinity of the reaction plane. It is, how- 
ever, possible that the physics backgrounds may be quite 
different for the same- and opposite-sign pairs, and the 
same- and opposite-sign difference may be dominated by 
physics backgrounds. For example, local charge conser- 
vation will naturally cause differences between the same- 
and opposite-sign pairs [29]. In fact, the results shown 
in Fig. |4] indicate that the centrality dependence of the 
asymmetry correlations is similar to the centrality de- 
pendence of the elliptic anisotropy. This is more clearly 
shown in Fig. [12] where the difference between the same- 
and opposite-sign results (scaled by -/Vp art ) is plotted as 
a function of the measured average elliptic anisotropy in 
each centrality bin. The dependence is roughly linear; 
the lines in Fig. 1121 show two linear fits, one with the in- 
tercept fixed at zero and the other with the intercept as a 
free parameter. If the charge separation is indeed a corre- 
lation background, then the approximate proportionality 
suggests that the charge-dependent correlation strength 
is insensitive to centrality. However, the apparent linear 
relationship does not necessarily mean that the charge 
separation must be an anisotropy related background. 
Because the CME and the average anisotropy are both 
functions of centrality, they can be indirectly related re- 
sulting in an apparent relationship between the charge 
separation and the average anisotropy. 

In order to gain further insights, one wants to fix the 
centrality hence the CME, and vary the event anisotropy. 
This can be achieved by the study in Fig. [7] of the asym- 
metry correlations as a function of the event- by- event 
elliptic anisotropy of the measured particles. Figure [7] 
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FIG. 12: The values of A(A 2 ) - A(A+A-) , scaled by iVpart, as 
a function of the measured average elliptic anisotropy s ) 
in Au+Au collisions. The centrality bin number is labeled by 
each data point, for 70-80% up to 8 for 0-5%. The error 
bars are statistical only. 



suggests, given a fixed range of centrality, that the bulk 
event structure may have a significant effect and the 
backgrounds for same- and opposite-sign pairs may in- 
deed be different. The results in Fig. [7] could be inter- 
preted as follows. The values of S(Al R ) decrease with 
increasing v% hs , while the values of S(Ay T) ) increase. The 
trends of S(A[ j - R ) could result from a relative abundance 
of back-to-back same-sign pairs in-plane rather than out- 
of-plane. The more abundant back-to-back pairs in-plane 
give a larger v% hs and reduce the LR asymmetry, thereby 
decreasing S(Al R ). Likewise, the S{A 2 }T) ) trends could 
result from a reduction in the back-to-back same-sign 
pairs out-of-plane rather than in-plane, which increases 
both the i>2 bs and <^(^ud)- The v% hs dependences in 
<5(A + A_)tjd and 6(A+A-)lr are significantly weaker. 
The trends seem to be opposite from those in 5(A 2 J - D ) 
and S(A 2 jR ). This may stem from the different nature 
of the correlations between opposite-sign pairs (small- 
angle) and same-sign pairs (back-to-back). These be- 
haviors of S(A 2 ) and 8(A+AJ) with v^° s may be in- 
line with suggestions that those charge correlations arise 
from cluster particle correlations overlaid with elliptic 
anisotropy (28l.[29j. 

Figure [13] (left panel) shows the difference between 
same- and opposite-sign correlations, A = A(A 2 ) — 
A(A + A-), as a function of the event- by- event v% hs in 
20-40% central Au+Au collisions. At large positive i>2 bs , 
A(A 2 ) > A(A + A^) is consistent with the CME. It is pos- 
sible that at significantly negative v% hs , the reconstructed 
EP may be orthogonal to, rather than aligned with, the 
real reaction plane so that UD and LR are flipped. As 
a result, the negative A would really be positive if cal- 
culated related to the true reaction plane. This would 
also be consistent with the CME. On the other hand, for 
events with modest negative V2 bs > —0.1, it is found by 
the sub-event method that the EP resolution is relatively 



well defined (see Fig. [2"5l in Appendix IB 6[) . However, in 
the region —0.1 < V2 bs ~ 0, the values of A are neg- 
ative. This suggests that the CME, which should give 
A(A 2 ) > A(A + A-) , cannot be entirely responsible for 
the present observations. 

The STAR data from RHIC Run VII have also been an- 
alyzed using the first harmonic event plane reconstructed 
by the ZDC-SMD [H H The ZDC-SMD event 

plane resolutions can be found in Appendix lB 61 The cor- 
responding same- and opposite-sign correlation difference 
as a function of the event- by- event elliptic anisotropy rel- 
ative to the ZDC-SMD event plane is shown in Fig. [TS] 
(middle panel). This result agrees with that obtained 
using the second-harmonic TPC event plane. 

Several authors [H, [2i| H3| have suggested that there 
could be a physics background proportional to the event- 
by-event ^2 bs due to the net effect of particle intrinsic 
correlations and production elliptic anisotropy. Thus, the 
values of A(«2 bs ) were fit with a straight line versus v% hs . 
The result is A(vf s ) = (1.3±1.4(stat) j^;°(syst)) x 10~ 5 + 
(3.2 ± 0.2(stat)tg;|(syst)) x lO^f 8 for 20-40% central 
Au+Au collisions. 

If the linear dependence of the data is entirely due 
to background correlations of the type suggested in 
Ref. [H, [H, [H and the CME does not contribute ap- 
preciably to the final state i>2 bs , then the intercept is the 
most sensitive to possible CME effects, and the slope 
is a measure of the background correlation strength. 
In this case, multiplicity asymmetries of particles in 
phase space which are less elliptically distributed (near 
«2 bs = 0) [41], with respect to an event plane recon- 
structed elsewhere in phase space, will be more sensi- 
tive to possible CME effects. Shown in Fig. [13] (right 
panel) is the charge separation A = (A 2 ) — A(A + A_) 
in events with |v2 bs | < 0.04 as a function of the cen- 
trality. The values of A are multiplied by iV par t for 
clarity. The charge separation seems to be consistent 
with zero within the present statistical precision, sug- 
gesting no substantial charge separation in those events 
with approximately zero ellipticity of the measured par- 
ticles. For the mid-central 20-40% events, it was found 
that A(|v° bs | < 0.04) = (1.9±1.9(stat)±|j(syst)) x 10~ 5 . 
For comparison, the charge separation averaged over all 
events is A = (1.6 ± O.l(stat)) x 10~ 4 . 

Table [I] lists the Run IV TPC data of the charge sep- 
aration parameter A, the charge separation parameter 
within the range of |v£ bs | < 0.04, and the linear fit in- 
tercept and slope as a function of centrality. The first 
quoted error is statistical and the second quoted error 
is systematic. The details of the systematic uncertainty 
study are given in Appendix [Cl 

To gain more insight, the charge asymmetry variances 
and covariances using UD— LR differences with respect 
to a random azimuthal plane as a function of the event 
elliptic anisotropy relative to the randomly chosen plane 
were also studied. The dependences of UD— LR versus 
«2 bs are shown in the left panel of Fig. [Inland are qualita- 
tively similar to those with respect to the reconstructed 
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FIG. 13: (Color online) A = A(A 2 ) - A(A+A-} as a function of v% bB , the event-by-event elliptical anisotropy of particle 
distributions relative to the second-harmonic event plane reconstructed from TPC tracks (left panel) and the first harmonic 
event plane reconstructed from the ZDC-SMD neutron signals (middle panel), in 20-40% central Au+Au collisions. Right 
panel: Average A for events with \v% s | < 0.04 relative to the TPC event plane as a function of centrality. The error bars are 
statistical only. 



TABLE I: The charge separation parameter A = A{A 2 } — A{A+A-) for all events and for events with the observed second 
harmonic parameter \v2 hB \ < 0.04, and the linear fit intercept and slope to A(i>2 bs ), as a function of centrality (.Apart is the 
corresponding number of participants). The first error is statistical and the second asymmetric error is systematic. 



Centrality 


Apart 


A 


A(\vf s 


< 0.04) 


intercept 


slope 






80-70% 


15.7 


(3.0 ± 1-8±|;|) x 10" 4 


(-3.0 ±4.1 


+9.4 
-2.4 


) x 10" 4 


(0.2 ±2.0±^) x 10" 4 


(1.4±0.2±g;«) 


x 10" 


2 


70-60% 


28.8 


(4.6 ± O.8+10) x 10~ 4 


(3.3±1.7t° 2 ; 4 ) 


x HT 4 


(2.4±0.9i^) x 10~ 4 


(8.6 ±0.915;?) 


x 10" 


3 


60-50% 


49.3 


(3.8 ±0.51^) x 10" 4 


(0.4±0.8j 


1.9\ 
2.5/ 


x HT 4 


(1.1 ±0.5t?;?) x 10_4 


(6.9±0.5±?;g) 


x 10" 


-3 


50-40% 


78.3 


(2.7±0.3j£?) x 10" 4 


(0.7±0.5l 


l.lN 

1.6) 


x 10" 4 


(0.4 ± 0.3±i o) x 10" 4 


(5.1±0.4i° ^) 


x 10" 


3 


40-30% 


117.1 


(2.0±Q.2ig;?) x 10" 4 


(2.8 ± 3.1 j 


7.2 x 
6.l) 


x 10 -5 


(1.5±2.2l[-j) x 10" 5 


(3.8 ± 0.3±g;J) 


x 10" 


3 


30-20% 


167.6 


(1.2±0.lj£?) x 10" 4 


(1.1 ± 2.0j 


3.1% 
3.6/ 


x 10~ 5 


(0.8±1.5li;g) x 10~ 5 


(2.6±0.2±° ; 4 ) 


x 10" 


-3 


20-10% 


234.2 


(5.8 ± 0.8^) x 10" 5 


(2.1 ± 1.3 j 


1.0 -I 
3.2^ 


x 10~ 5 


(1.4 ± 1-0±|;|) x 10 -5 


(1.4±0.2±° ; 4 ) 


x 10" 


-3 


10-5% 


298.6 


(2.6 ± 0.9ti;?) x 10" 5 


(-0.6 ±1.2 


+3.2 
-0.9 


) x 10 -5 


(0.5 ± l.Otj^) x 10" 5 


(1.1 ± 0.2±g;J) 


x 10" 


3 


5-0% 


350.6 


(-0.3 ± O.Slo.g) x 10~ 5 


(-1.7 ± 1.0 


+3.2 
-0.4 


) x 10" 5 


(-1.3±0.8±^) x 10~ 5 


(8.5±2.2±?; 4 ) 


x 10" 


4 



EP. The A(A 2 ) and A(A + A_) from the random plane 
cross at vf 311 ~ with approximately zero intercept, while 
those from the recontructed EP cross at a positive in- 
tercept as aforementioned. The differences between the 
reconstructed EP and random plane results are shown in 
the middle panel of Fig. Q31 The differences are sensitive 
to EP-dependent correlations, and appear to be indepen- 
dent of charge signs. The difference between the same- 
and opposite-sign charge UD— LR fluctuations, shown by 
the open triangles in the right panel of Fig. [TJ1 are also 
linearly dependent on the observed i>2 bs relative to the 
random plane. The dashed line is a linear fit to the open 
data points. For comparison the data with the recon- 
structed EP and the corresponding linear fit from the 
left panel of Fig. [13] are superimposed as the solid points 
and solid line, respectively. The linear dependences of 
the charge separation A(«2 bs ) are equal between the re- 
constructed EP and the random plane within the sta- 
tistical uncertainties. This confirms that the observed 
charge separation is indeed correlated with the observed 
final-state event shape. 

It is worthwhile to note that, with respect to the ran- 



dom plane, the average (i>2 bs ) is zero (the event proba- 
bility distribution is symmetric about v% hs = 0), and the 
event-integrated charge separation is zero as expected. 
With respect to the reconstructed EP, the center-of- 
gravity of the event probability distribution is no longer 
at zero but at positive (v% hs ) (i.e. the observed average 
elliptic anisotropy) and the event-integrated charge sep- 
aration is finite and positive as reported in the present 
work. 



Further comparisons of oriented and random event 
plane results will be interesting. Additional studies to 
disentangle effects from the final-state event shape and 
the magnetic field have been proposed by using collisions 
of uranium nuclei [42j which have a large intrinsic defor- 
mation. These studies will provide complementary exper- 
imental data for understanding the origin of the observed 
charge separation. 
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FIG. 14: (Color online) Left panel: A(A 2 ) and A{A+A-} as a function of with a random EP in 20-40% Au+Au collisions. 
The lines depict the results with reconstructed EP from Fig.[7]lower right panel for comparison. Middle panel: The differences in 
A(A 2 } and A(A+A-} between the reconstructed EP and random EP results, respectively. Right panel: A = A{A 2 } — A{A+A-) 
as a function of i>2 bs with a random EP (open triangles). The solid triangles are results with reconstructed EP from Fig. [T3] 
left panel for comparison. Error bars are statistical. 



VI. SUMMARY 

Correlations of positive (A + ) and negative (A J) charge 
multiplicity asymmetries with respect to the event plane 
and the plane perpendicular to the event plane have been 
measured. The asymmetries are measured in one half of 
the STAR TPC (one unit in pseudo-rapidity) while the 
event plane is reconstructed in the other half of the TPC. 
The dynamical variance 6(A±) and covariance 5(A + A-) 
are obtained by subtracting the contributions from the 
statistical fluctuations and detector effects, and were pre- 
sented in d+Au and Au+Au collisions. 

The charge asymmetry dynamical variances, S(A 2 j T) ) 
and S{Aj jR ), are positive in d+Au and peripheral Au+Au 
collisions and become negative in medium-central to cen- 
tral collisions (cf. Fig. [3]). The positive dynamical vari- 
ances in ei+Au and peripheral Au+Au collisions indi- 
cate that the dynamical physics processes broaden the 
charge asymmetry distributions. The same-sign pairs 
are preferentially emitted in the same direction. The 
negative dynamical variances in medium-central to cen- 
tral collisions indicate that the dynamical physics pro- 
cesses narrow the charge asymmetry distributions. The 
same-sign pairs are now preferentially emitted back-to- 
back. The charge asymmetry covariances, 5(A + A-)\je, 
and 5(j4 + A_)lr, are positive and large, indicating that 
significant correlations are present between positive and 
negative particles in the events (cf. Fig. [3]). The posi- 
tive and large covariances indicate that the opposite-sign 
pairs are preferentially emitted in the same direction, and 
the aligned emission is largely independent of the reac- 
tion plane. 

The charge asymmetry dynamical variances, 6(A 2 j- D ) 
and 5(A 2 jB ), decrease with increasing pr (cf. Fig. 
They are positive at low px, indicating that the \ow-px 
same-sign pairs are preferentially emitted in the same 
direction. They are preferentially back-to-back at large 
Pt in central collisions. The charge asymmetry covari- 
ances are approximately independent of pt at low pt, 



and rapidly increase with px- The aligned same-side 
emission of opposite-charge pairs is significantly stronger 
at large pr- 

The dynamical charge asymmetry variance, S(A^ ! - D ), is 
larger than S(Af j j t ). The charge asymmetry covariance, 
8(A + A-)ij-d, is larger than 5(/1_|_^4_)lr in non-peripheral 
Au+Au collisions. The differences A(A 2 ) = 5(A 2 ]T) ) - 
S(Al n ) and A(A+A_) = S(A + A-) UU - 6(A+A-) LR in- 
crease with centrality in peripheral collisions, reaching 
a maximum in medium central collisions, and then de- 
creases towards more central collisions (cf. Fig. 0}. The 
differences, A(A 2 ) and A(A + A-) , in a given centrality 
bin increase with px (cf. Fig. The centrality and px 
dependences are qualitatively similar between the same- 
and opposite-sign UD— LR correlations. 

To gain more insight into the contributing physics 
mechanisms, two aspects of the charge asymmetry cor- 
relations were investigated. First, the charge asym- 
metry correlations as a function of the event- by- event 
elliptic anisotropy of the particles in the asymmetry 
measurements, v% hs , was studied (cf. Fig. [7|). The 
same-sign UD— LR correlation, A(A 2 ), increases with 
increasing v% hs , and the opposite-sign UD— LR correla- 
tion A(A + A^) decreases slightly. They cross at the 
same positive value at U2 bs ~ 0. Also studied was 
the charge asymmetry correlations as a function of the 
size of the azimuthal region for the asymmetry measure- 
ments (cf. Fig. [8|). It was found that the UD— LR dif- 
ferences of both the variance and covariance, A(A 2 ^^) 
and A(v4 +j 4_)a0 w , increase with decreasing wedge size. 
However, the difference between the two diminishes with 
decreasing wedge size. 

The present measurements were motivated by the Chi- 
ral Magnetic Effect, which could yield a charge separation 
of final state hadrons across the reaction plane. Previ- 
ous STAR measurements of the three-particle correlators 
were consistent with charge separation by the CME to- 
gether with medium interactions. However, the possi- 
ble physics backgrounds for the measured correlator re- 
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suits were not fully explored. Similarly, the CME charge 
separation would yield a positive same-sign correlation 
A(A 2 ) and a negative, or zero in the case of medium 
interactions, opposite-sign correlation A(A + A-). The 
present measurements show A = A(A 2 ) — A(A + A_) > 0, 
an indication of charge separation, although both A(A 2 ) 
and A(A + A-) arc positive. For the mid-central 20-40% 
events, it was observed that A = (1.6±0.1(stat)) x 10~ 4 . 

It is possible that the physics backgrounds for same- 
sign and opposite-sign correlations are equal and fall in- 
between the present UD— LR results. In such a case, only 
the difference A = A(A 2 ) — A(A + A-) is sensitive to the 
CME and their average behavior is due to other physi- 
cal mechanisms. The measured values of A were largest 
when whole hemispheres were used in the asymmetry 
measurements and diminishes when smaller wedge sizes 
were used (cf. Fig. [TU]). This suggests that the charge 
separation across the event plane is mainly in-plane and 
same-sign pairs are roughly back-to-back in the upper or 
lower hemisphere. 

It is also possible that the physics backgrounds are 
different for the same- and opposite-sign UD— LR corre- 
lations. In fact, the charge separation, A = A(A 2 } — 
A(A + A-), is roughly proportional to the average ellip- 
tic anisotropy in each centrality bin (cf. Fig. IT2]> . It 
is further found that the values of A are approximately 
proportional to the event-by-event elliptic anisotropy of 
the particles measured in the charge multiplicity asym- 
metries: A(vf s ) = (1.3 ± 1.4(stat)^;°(syst)) x 10" 5 + 
(3.2 ± 0.2(stat)±£|(syst)) x lCT 3 ^ 8 This proportional- 
ity was also observed using the first-order harmonic event 
plane reconstructed using the ZDC (cf. Fig. IT3]) . This 
suggests that the physics backgrounds may be a net ef- 
fect of particle intrinsic correlations and production ellip- 
tic anisotropy. The intrinsic particle correlations include 
momentum conservation and local charge conservation, 
although the charge conservation effects are only present 
for opposite-sign charges. In events of nearly zero par- 
ticle ellipticity (|f2 bs | < 0.04) where such a background 
may be absent, the charge separation effect was observed 



to be A(\vf s \ < 0.04) = (1.9±1.9(stat)+^|(syst)) x 10 
for the 20-40% centrality. 
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Appendix A: Connections to Correlators 

The differences (A| UD ) - (-A±. L r) and (A + A_) UD - 
(A|_A_)lr are related to the previously reported three- 
particle azimuth correlators fLU [l5| , but with important 
differences. By using a Fourier series of a step-function 
in <fi — ipsp (the particle azimuth relative to the event 
plane) , the present charge asymmetry observables can be 
expressed as 

a 4 ^ sin(2n + l)(4>±, l - V'ep) 



.4 



±,LR 



ttN± 
4 



j=l n=0 
N± oo 



2n+ 1 



y^ y^ cos(2n + l)(</>±,» - V'ep) ^^ 



2n+ 1 



1=1 n=0 

The difference in asymmetry correlations is 



2 / N a ,Np oo 



N a Np .4-' ^ n (2n + l)(2m + l) 

4 ^ i,j=0 n.rn—0 \ / \ / 

where a and f3 stand for or '— 5 particles. This is similar in form to the azimuth correlator observable (co$((j) a 4- 
4>f3 - 2^ RP )) « (cos(0 a + fa - 20 c ))/« 2 ,c of Ref. 0, HjJ: 

(cos(4> a + 4>/3- 2V>rp)) = ( N , N _ s ^2 cos(0 Q +<t> p - 2tp RP ) \ , (A3) 



I 

where ^>rp is the reaction plane angle and 8 a p is the contain all possible (including mixed) harmonic terms, 
Kroncckcr delta. The asymmetry correlation differences while the correlator observables contain only one of the 
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infinite number of terms. While the correlators are in 
terms of azimuthal angle relative to ipBP , the asymmetry 
correlation observables are in terms of azimuthal angle 
relative to V>ep, and are therefore affected by EP resolu- 
tion. 

The present charge asymmetry correlations and the az- 
imuthal correlators are related but differ as shown in 
Eq. (|A2I) and Eq. (|A3[) . To gain more insight into the 
relationships and differences, the upper panel of Fig. 1151 
shows the comparison between the asymmetry correla- 
tion differences (LR — UD) and the azimuthal correlator 
(cos(0 Q + <j)p — 2^rp)). The charge asymmetry correla- 
tion differences LR — UD in Fig. [15] have already been 
divided by the event-plane resolution. Those correlators 
are calculated using the a and ft particles from the same 
r\ region as were used for the charge asymmetry measure- 
ment with particle c from the other r\ region used for EP 
construction. The correlator results shown in Fig 1151 are 
not identical to those published previously 0, [l5[ where 
the three particles a, ft, and c are from the entire TPC 
acceptance of —1 < r\ < 1. The present correlator mea- 
surements, which are from a smaller range in tk hence a 
smaller range in Arj, are larger than those in [lj, Ell be- 
cause the measured correlators decrease with increasing 
range of At? (HQjI. 

To compare directly to the azimuthal correlators in 
Refs. [lj, Il5} . the present charge asymmetries were 
formed using the entire TPC acceptance. To avoid the 
self-correlation, the event was divided randomly into two 
sub-events. One sub-event was used to reconstruct the 
event plane and the other was used to calculate the 
charge asymmetries. The results are shown in the lower 
panel of Fig. 1151 a nd are compared to the azimuthal cor- 
relators from |l4, HE)]. The azimuthal correlators have 
also been calculated using two particles, a and b, from 
the sub-event used for the asymmetry measurements and 
the third particle c from the other sub-event used for 
EP reconstruction. The azimuthal correlators (shown 
as open points in Fig. IT5"j) are consistent with those of 
Ref. QES. 

As shown in Eqs. (|A2[) and (|A3I) . the correlators con- 
tain an infinite number of harmonic terms in the asym- 
metry correlation difference observables. The present 
values of (Al R ) — (A^ D ) are comparable to the same- 
sign correlator (cos((/> Q + 4>p — 2?/>rp)). This suggests 
that the high-order harmonic terms in A(A 2 ) may be 
small. This in turn suggests that the event-plane res- 
olution correction by linear extrapolation may be suffi- 
ciently accurate for the measurements of (^4lr) ~ (^ud)- 
On the other hand, the present values of (A + A_)lr — 
(A + A_)ud is significantly different from the opposite- 
sign correlator. It is important to note that the differ- 
ence (^4 +j 4_)lr — (A +j 4_)ud is significantly negative, 
while the correlator magnitude is close to zero. This 
suggests that the high-order terms in Eq. (|A2I) are im- 
portant. Because of these terms, the event-plane reso- 
lution correction via the linear extrapolation applied on 
(A+A_)lr— (^4+^4-)ud in Fig.[T5]is likely invalid. How- 
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FIG. 15: (Color online) The correlation differences, (v4lr) — 
(^ud) an d (^4+^4-}lr — (A+A_)ud, scaled by the number 
of participants (7r/4) 2 A'p ar t. Also shown are (cos(cj> a + 4>p — 
2(f> c ))/v2, c ~ {cos(4> a + 4>p — 2-i/>Rp)} of same-sign and opposite- 
sign particle pairs (a and ft) calculated from the same 77 region 
as the charge asymmetry correlations with particle c from the 
other rj region used for EP construction. The error bars are 
statistical only. 



ever, the imperfect event-plane resolution can only reduce 
the magnitude of (A + A_)ud — (^4+^4-)lr that is shown 
in Fig. H Thus, the true (A + A-) LR - (A + A-) UD with 
respect to the real reaction plane may be even more neg- 
ative than that shown. 



Appendix B: Analysis Details 

1. Detector Efficiency 

Figure [TH] (upper panel) shows examples of the (jy 
distributions of the positively charged particle multiplic- 
ity within the 77 > and 77 < regions separately for 
30-40% central Au+Au collisions. The negative particle 
results are similar. The results from other centralities are 
also similar. In Fig.[T51 the magnetic held polarities have 
been summed and the particles are integrated over the 
range 0.15 < pt < 2 GeV/c. The regular-pattern of the 
TPC sector boundaries is clearly seen for the both posi- 
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tive and negative rj particles. The rj < particles have an 
additional inefficiency in the region of ll7r/6 < <fi < 2tt, 
due to an inefficiency in the electronics for two of the 
sectors in the east side of the TPC. This inefficiency per- 
sisted over the whole period over which the present data 
was collected but showed no significant variation over 
time within this period. 




The efficiency corrections were obtained in the follow- 
ing way. The multiplicity distributions in </3, such as 
the ones in Fig. [T5] (upper panel), were normalized to 
average unity. These normalized distributions are re- 
ferred to as the acceptance x efficiency, e(4>). These 
e(4>) distributions are separated according to the mag- 
netic field polarities, particle charge signs, and centrality 
bins. They are further separated for positive and nega- 
tive r] (approximately corresponding to the two halves of 
the TPC depending on the primary vertex position), and 
for the following pt bins: 0.15-0.5, 0.5-1.0, 1.0-1.5, and 
1.5-2.0 GeV/c. The detector acceptance x inefficiency 
correction factor is taken as 1/e (</>). The correction fac- 
tor l/e((j>) here can be larger or smaller than unity. 

The overall single particle tracking efficiency is not 
corrected for during the event-plane construction. The 
present results were checked when using an event-plane 
reconstruction for which the particle multiplicities were 
corrected by the centrality and p-r-dependent tracking 
efficiency. This did not affect the present results. 




FIG. 16: (Color online) Upper panel: Positively charged par- 
ticle multiplicity distributions versus the azimuthal angle, <j>, 
in 30-40% central Au+Au collisions. The data are shown sep- 
arately for r) > and r\ < 0. The magnetic polarities of the 
STAR magnet have been summed and the pr is integrated 
over the range 0.15 < pr < 2 GeV/c. The inverse of these 
distributions, properly normalized, were used to correct for 
the track efficiency versus (j>. Lower panel: Reconstructed 
event-plane azimuthal angle distributions in 30-40% central 
Au+Au collisions. Particles within 0.15 < pr < 2 GeV/c 
from rj < and r\ > were used separately to reconstruct the 
EP. The error bars are statistical only. 

The single particle (^-dependent inefficiencies, primar- 
ily due to the TPC sector boundaries, were corrected for 
in reconstructing the EP and in calculating the multi- 
plicity asymmetries. Figure [16] (lower panel) shows the 
event-plane V'ep distributions reconstructed separately 
from the r\ > and rj < particles in 30-40% central 
Au+Au collisions. The distributions are uniform. Fit- 
ting the 1] > and r\ < event-plane distributions with 
a constant value resulted in values of x 2 /ndf of 39/39 and 
48/39, respectively. 



2. Self-Correlation 

Four asymmetries were calculated separately, (i) using 
particles within < 77 < 1 with the EP reconstructed 
from — 1 < rj < 0, (ii) using particles within — 1 < rj < 
with the EP reconstructed from < rj < 1, (hi) us- 
ing particles within — 1 < r\ < with the EP recon- 
structed from — 1 < rj < 0, and (iv) using particles within 
< rj < 1 with the EP reconstructed from < 77 < 1. 
The UD charge asymmetry covariance relative to that of 
LR, (A + A_) UD /(yl + A_) LR , is shown in Fig. [17] for all 
four cases. Differences are observed between those us- 
ing the same set of particles for the EP construction and 
asymmetry calculation (cases (iii) and (iv)) and those us- 
ing different sets of particles (cases (i) and (ii)). The dif- 
ferences are more significant in peripheral collisions. The 
reason for this difference is as follows. The reconstructed 
EP divides the multiplicity of the event into two roughly 
evenly populated halves. Therefore, the positive charged 
particle asymmetry and that for negative particles cal- 
culated using the same set of particles used for the EP 
reconstruction are anticorrelated between UD. This does 
not affect those of LR. This results in a relatively smaller 
(A + A^)uv than (A + A_)lr for cases (iii) and (iv). To 
avoid this self-correlation, data from cases (i) and (ii) 
were used, where the particles used for the EP construc- 
tion and for the asymmetry calculation are different. On 
the other hand, no apparent self-correlation is observed 
in the variances (A± UD ) and {A 2 ± LR ). This is because, 
in the cases of (iii) and (iv), only half of the particles 
used in the EP reconstruction are used for asymmetry 
calculation. However, cases (iii) and (iv) are not used in 
the analysis of the variances (A 2 ± UD ) and (A\ LR ) to be 
consistent with that of the covariances (/!+/!_ )ud and 
(A+A_) LR . 
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FIG. 17: (Color online) The relative charge asymmetry corre- 
lations, ( j 4+ j 4_)ud/(^1+^4-}lR; as a function of the number 
of participants, N pax t, for four combinations of r\ ranges used 
for EP reconstruction and asymmetry calculation. The error 
bars are statistical only. 



3. ^-Efficiency Correction 

The TPC sector boundaries could create "dynami- 
cal" event-by-event correlations even though the EP an- 
gles are uniform. Specifically, the east TPC sector with 
systematically lower efficiency introduces a "dynamical" 
event- by- event correlation (see below). To reduce this 
sort of fake dynamical correlation, the sector boundary 
effects were corrected for by the (/(-dependent correction 
l/e(cj>). This correction is applied on average over the 
entire dataset. In principle, any time variation in e((/>) 
would introduce "dynamical" event-by-event correlations 
if a single set of correction factors is applied. However, it 
was found that e(<p) is stable over the entire run period 
of the present data. 

Figure [TS] shows the charge asymmetry correlations 
ud) an d (^4+^4-)ud before and after the 
^-efficiency corrections. The asymmetry correlations are 
multiplied by the number of participants, N paYt , to better 
show the magnitude. For clarity, only the r\ > region is 
shown for (A 2 ^ UD ) and only the 77 < region is shown for 
{A 2 _ UD ). The values of (A 2 , UD ) and (A 2 _ UD ) are similar 
in the same rj > or 77 < region. The LR asymmetry 
correlations are similar to those of UD. The corrections 
for 77 < are larger than that for 77 > 0. This is due 
to the greater non-uniformity from the TPC electronics 
inefficiency at ll7r/6 < < 27r in the 77 < region (see 
Fig.HU). ' 

As seen in Fig. IT51 (right panel), before the (^-efficiency 
correction, the values of (A + A-) for the 77 < region 
are significantly larger than those in the 77 > region in 
central collisions. After the ^-efficiency correction, the 
values of (A + A_)tjd from the 77 < and 77 > regions 
are consistent, as will be shown in Fig. [20] (right panel). 
As for (A 2 ), the corrections for (A + A-) for 77 < are 
relatively large. In fact, the inefficient sector boundaries 



seem to have similar effects on the absolute magnitudes 
of (A 2 ) and (A + A-). Relatively, they are more signifi- 
cant on (A + A-) than on (A 2 ). The values of (v4 +j 4_)lr 
are not shown in Fig. [TS] (right panel) but they are sim- 
ilar to (A + A_)tjd- The curves shown in Fig. [IS] will be 
discussed in Appendix IB 41 



4. Statistical Fluctuation and Detector Effects 

The variances are non-zero even when no dynamical 
fluctuations are present. This is due to the trivial effects 
of statistical fluctuations of the multiplicity. If one takes 
A^±,u = (N± t v) + 6N± )XJ and A±, D = (JV ±)D ) + <5A±, D 
(where N± v collectively stands for N +y \j and N- t \j, and 
N±j3 for N + p and N- t u), and assume the fluctuations 
are Poissonian, the statistical fluctuations can be ex- 
panded in 5N±/(N±) and approximated by 



(^±,UD,stat) 




5N±y - 5N±p 
(SN ±tV + 6N±p)/(N±] 



(B 



A similar formula exists for (A± LRstat ). The multiplici- 
ties in Eq. (|B1[) are the measured multiplicities, not the 
efficiency-corrected ones. The average efficiency correc- 
tions cancel in the present charge asymmetries and do 
not contribute to the statistical fluctuations. The result 
by Eq. (|B1[) for (A^_ LR stat ), as an example, is shown in 
Fig. [19] (left panel) as the dotted curve. This approxi- 
mation using Eq. (|B1[) will be referred to as the "1/iV" 
approximation. 

In the data analysis, one can obtain the contributions 
from statistical fluctuations by assigning each particle 
randomly into the up or down hemisphere and into the 
left or right hemisphere (referred to as "50-50" method), 
and then calculating the asymmetry correlation without 
additional corrections by l/e(0). This would correspond 
to the case of purely statistical fluctuations with a per- 
fect detector. Figure Q1J] (left panel) shows the values of 
{A 2 + LR stat ) in the region 77 < obtained by this "50-50" 
method as the dashed curve. As seen from the figure, 
the l/N approximation of Eq. (|B1I) underestimates the 
statistical fluctuations. To quantify the magnitude of the 
underestimation, Fig. Q1J] (middle panel) shows the ratio 
of the dashed to the dotted curves from the left panel. 
Also shown is the corresponding ratio from the 77 > re- 
gion. As expected, the statistical fluctuations relative to 
the l/N approximation are the same between the 77 > 
and 77 < regions. The underestimation is the most se- 
vere in peripheral collisions, and is of the order of a few 
percent. This is presumably due to the non- Poissonian 
nature of the multiplicity distributions in the peripheral 
data and the approximations made in Eq. (IB1|) . In cen- 
tral collisions, both the Poissonian multiplicity distribu- 
tion assumption and the large N approximation should 
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FIG. 18: (Color online) The asymmetry correlations, (y4+ iUD ) (left panel), (At, ud) (middle panel), and (^4+^4-)ud (right 
panel), multiplied by the number of participants iVpart, before and after the corrections for the ^-dependent acceptance x 
efficiency. The colored curves are the corresponding statistical fluctuations and detector effects. The results are separated for 
■q > and r\ < 0. Only 77 > is shown for (A+ UD ) and r? < for (At. UD ) for clarity. The error bars are statistical only. 
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FIG. 19: (Color online) Left panel: Statistical fluctuation and detector effects in the charge asymmetry variance (multiplied 
by the number of participants iVpart) from the r\ < region. The dotted curve shows the 1/N approximation using Eq. flBll l. 
the dashed curve shows the statistical fluctuations (AL LR., s tat) obtained by the "50-50" method (see the text). The solid curve 
connecting the solid points shows the net effect of the statistical fluctuations and detector non-uniformities, (At. lr stat+dct)i 
obtained by the adding-7r method, and the crosses shows the same but obtained from the "scramble" method. See the 
text for details. Middle panel: The statistical fluctuation effects relative to the 1/N approximation (thin curves) and the 
effects due to the imperfect detector (thick curves). The dashed curves are for the region 77 > and the solid curves are 
for r\ < 0. Right panel: The ratios of the statistical fluctuation and detector effects (^4+.LR,stat+d c t)/(^-.LR,stat+d c t) an d 
(At_ uD,stat+dct)/( j 4+,LR,stat+dct)i separately for the 77 > and r\ < regions. The (^-acceptance correction was applied for the 
results in this figure. The error bars are statistical only. 



be valid. As can be seen, the statistical fluctuations in 
central collisions can be well described by Eq. (|B1|) . 

The STAR TPC is not precisely uniform in azimuth. 
There are sector boundaries and different detection effi- 
ciencies due to variations in the electronics performance. 
As mentioned earlier, these non-uniformities are cor- 
rected for on average by the ^-dependent correction fac- 
tor l/e(4>) separated by particle charge signs, magnetic 
field polarities, and positive and negative rj regions, cen- 
trality bins, and pt bins. However, the event-by-event 
fluctuations in the efficiencies cannot be corrected. These 
fluctuations in detector performance introduce a dynam- 
ical effect, which are referred to as "detector effects" in 
this paper. These detector effects were quantified in the 
following way. For each particle, one either adds ir to 



the measured azimuthal angle <f> or does nothing. The 
probability to add n is determined by the relative effi- 
ciency at 4> and + namely, e((f) + n) / {e{<j>) + e((f) + n)) . 
In this way, the physics correlation amongst particles is 
destroyed, but the detector non-uniformities and elliptic 
flow correlations of the particles are preserved. The av- 
erage (^-dependent efficiency is then corrected for each 
particle depending on its new (f> values, and the charge 
asymmetry correlation were then calculated. The resul- 
tant asymmetry correlation from the region 77 < for 
positive particles, (A 2 ^ LR stat+dct ), is shown in Fig. [TTfl 
(left panel) by the solid curve connecting the dots. The 
text "stat+det" is used to stand for the net effect of the 
statistical fluctuations and the detector non-uniformities. 
The solid curves in Fig. [T5] (middle panel) show the ratio 
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of ( A +,LR, sta t+dct)/(^+XR, s tat)- This ratio measures the 
magnitude of the detector effects relative to the pure sta- 
tistical fluctuations. As seen from the figure, the ratio is 
approximately unity in peripheral collisions and increases 
to a few percent above unity in central collisions. This 
is expected because the detector non-uniformity is the 
strongest in the most central collisions [43(. The same 
ratio for the region 77 > is also shown. 

An alternative approach was also pursued. For each 
particle, a </> angle was randomly generated according to 
the efficiency e(<p). The charge asymmetry correlations 
for this "scrambled" data were then calculated in the 
same way as was done for the real data. This should 
be equivalent to the adding- 7r method above, except that 
it also destroys the contribution from elliptic flow. The 
results obtained by this "scramble" method are shown by 
the crosses in Fig. [19] (left panel) and are consistent with 
those obtained by the "adding- 7r" method (solid circles). 

Eq. (|B1I) indicates that, if there are no dynamical fluc- 
tuations, the self-correlation variable is approximately 
equal to the inverse of the particle multiplicity. Fig- 
ure [T2] (right panel) shows the ratio of the positive to neg- 
ative charge asymmetry correlations from the adding-7r 
method. Indeed, the (A 2 +stat+dct ) and (41 stat+dot ) dif- 
fer slightly. The values of (A 2 + stat+dot ) are smaller than 
those of (A 2 _ stat+dct ) by about 2%. This difference is 
consistent with the slight (about 2%) excess of positive 
particles over negative particles [221 ] . 

The statistical fluctuation and the detector effects 
should not depend on the reaction plane orienta- 
tion and they should be equal between UD and 
LR. It has been verified that (A% >m . stat+dct ) = 
(Al 

,LR, stat+dct 

) and (A 

-,UD,stat+dct — ,LR,stat+det/ ' 

Figure [19] (right panel) shows as examples the ratios 

Of (^+,UD,stat+dct)/(^+,LR, S tat+dot) fr0m the adding-TT 

method for the 77 > and r\ < regions. They are con- 
sistent with unity. The same is observed for the negative 
particle asymmetry correlations. 



5. Consistency Checks 

As shown in Fig. 1181 the charge asymmetry variances 
and covariances without applying the ^-acceptance cor- 
rections differ from those with the corrections. This sec- 
tion checks the corresponding statistical fluctuation and 
detector effects with and without the ^-acceptance cor- 
rections. The adding-7r method was used to calculate 
the statistical fluctuation and detector effects. The re- 
sults are superimposed as the curves of the corresponding 
colors in Fig. [18] for the 0-efficiency corrected and uncor- 
rected asymmetry correlation data. As seen, the statisti- 
cal fluctuation and detector effects for (A + AJ) without 
the (/^efficiency corrections are no longer zero, especially 
for the rj < region. However, the differences in the data 
points with and without the ^-efficiency corrections ap- 
pear to be the same as those between the corresponding 



curves for the statistical fluctuation and detector effects. 
This indicates that the differences observed in the data 
from different 77 regions are due to the statistical fluctu- 
ation and detector effects. This is demonstrated more 
quantitatively by the dynamical asymmetry correlations 
after subtracting the statistical fluctuation and detector 
effects. 

Figure [50] (left and middle panels) shows the in- 
efficiency corrected charge asymmetry variances (A 2 ^ UD ) 
and {A 2 _ UD ) from the 77 > and 7/ < regions 
separately. Superimposed are the combined statisti- 
cal fluctuation and detector effects (A 2 , UD stat+dot ) and 

(^-,UD,stat+dot)- The ( A Ir) are similar to the (A^). 
Figure [5D] (right panel) shows the ^-efficiency corrected 
covariances, (A + A^)\jd. The statistical fluctuation and 
detector effects in the covariances are zero, as indicated 
by the superimposed curves (obtained from the adding- ir 
method). The (j4+A_)lr is similar to the (A + A_)ud- 
In Fig. 1201 the minimum-bias d+ Au data are also shown. 
The analysis procedure for d+Au is the same as for 
Au+Au collisions. The d+Au data follow the trend of 
the peripheral Au+Au data. 

The asymmetry correlation results from 77 < and rj > 
should ideally be equal in Au+Au collisions because 
of the collision symmetry about mid-rapidity. This is 
approximately the case for the covariances, and both are 
consistent with zero. It is, however, not true for the 
variances. The corrected variances in rj < are larger 
than the corresponding ones in rj > by 1-3%. This is 
partially due to the fact that the variance magnitudes are 
dominated by the statistical fluctuations of multiplicities, 
which are approximately inversely proportional to the 
average multiplicity (see below). The 10% inefficiency in 
the ll7r/6 < </> < 27r region for 77 < results in a 2% 
larger magnitude for the variances compared to rj > 0. 
This is indeed shown by the two differing curves which 
are the corresponding statistical and detector effects. 

The above effects are cancelled in the difference be- 
tween UD and LR asymmetry correlations, which is the 
main observable sensitive to the local parity violation. 
This is because the ^-dependent acceptance corrections 
and the effect of statistical fluctuation of the multiplici- 
ties are identical for UD and LR asymmetries. 

The statistical fluctuation and detector effects 
{A\ s tat+dot) ( as obtained by the adding tt method) were 
subtracted from (A±) to obtain the dynamical asymme- 
try correlations: 8{A 2 ±) = (A 2 ± ) - (A^_ stat ). Figure [2T1 
(left and middle panels) show, respectively, the dynami- 
cal 5(A 2 [ _ UD ) and S(A 2 _ UD ) from the 77 > and 77 < re- 
gions, by taking the differences between the data points 
and the corresponding curves in Fig. [18] (left and mid- 
dle panels). Similarly, Fig. [5T] (right panel) shows the 
dynamical S(A+A-)ud. The corresponding statistical 
fluctuation and detector effects were subtracted from the 
charge asymmetry variances and covariances, calculated 
both with and without applying the (^-acceptance cor- 
rections. The resultant dynamical results are consistent 
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FIG. 20: (Color online) The asymmetry correlations, (A\ ud) (left panel), (A 2 _ UD ) (middle panel), and (A+A-)ud (right 
panel), multiplied by the number of participants N pax t, separately for r\ > and rj < 0. The star and triangle depict the results 
from d+Au collisions. The net effects of the statistical fluctuations and detector non-uniformities are shown as the curves. 
Error bars are statistical only. 
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FIG. 21: (Color online) The statistical and detector effect subtracted asymmetry correlations, 5(A\ LR ) (left panel), 5{A 2 _ LR ) 
(middle panel), and 8{A+A-)i,r (right panel), multiplied by the number of participants JVp art , before and after the corrections 
for the ^-dependent acceptance x efficiency. This figure corresponds to Fig. 1181 The error bars are statistical only. 



with and without the ^-acceptance corrections, as shown 
in Fig. [3TJ This indicates that no residual effects of the 
^-acceptance remains after the corrections for the sta- 
tistical fluctuations and detector effects. Nonetheless, 
the (^-acceptance corrections were applied in the present 
analysis. Unless specified, all results in this paper have 
the (^-acceptance corrections applied. 

A test was also performed where an acceptance hole 
was articifically created within a restricted <fi region by 
randomly discarding 50% of the particles in the <p region 
and at the same time reducing the acceptance x efficiency 
e(cj>) by a factor of 2. The analysis was repeated with 
the remaining particles to calculate the charge asymme- 
try variances and the statistical fluctuation and detector 
effects to the dynamical variances. The dynamical vari- 
ances were consistent with the results previously obtained 
without the artificial acceptance hole. This is also true 
for the charge asymmetry covariances. 

Figure |2"21 shows the dynamical 5(A 2 ± LR ) values from 
the rj > and rj < regions with the (^-acceptance cor- 
rections. This plot demonstrates that the results from 
the two rj regions are consistent, and the results are the 



same for positive and negative charge asymmetries. This 
is also true for UD asymmetries. 

Since the dynamical charge asymmetry variances from 
the rj < and r\ > regions are consistent, and 
those of the positive and negative charge asymmetry dy- 
namical variances are equal within the statistical un- 
certainties, as shown in Fig. [22j the average of the re- 
sults from the two rj regions is used. Also, the aver- 
age between the positive and negative charge asymme- 
try variances: S{A^ JT) ) = (8{A 2 + UD ) + 6 (At tjd))/ 2 an d 
S(At R ) = (S(A^ W ) + (5«lr))A is used. Since the 
covariances from the two rj regions are also the same (see 
Fig. [20] (right panel) and Fig. [21] (right panel)), the av- 
erage between the two regions is also used. The average 
variances and covariances are reported. 



6. Effect of the Event-Plane Resolution 

The constructed event-plane is not as same as the re- 
action plane. The inaccuracy, or event-plane resolution, 
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FIG. 22: (Color online) Dynamical charge asymmetry vari- 
ance after removing the statistical and detector effects, 

s ( a ±,lb) = <^±,lr> - (^±,LR, sta t+dot}> scale d by the num- 
ber of participants iVpart. The results are consistent between 
positive and negative 77 regions and between positive and neg- 
ative charges. Error bars are statistical only. 



is due to the finite multiplicity of particles used for the 
event-plane reconstruction. The observed differences be- 
tween the UD and LR charge asymmetry correlations are 
affected by the event-plane resolution. The magnitudes 
of the measured asymmetry correlation differences are re- 
duced from their true values (i.e. with respect to the real 
reaction plane) due to the finite event-plane resolution. 

The event-plane resolution can be calculated, approx- 
imately, by e EP = y/ (cos2(V>EP,r ? >o - V'EP,r ; <o)), where 
"0ep,?7>o and ?/>ep,tko are ^ ne reconstructed event plane 
azimuths from particles at rj > and r\ < 0, respectively. 
Note, this event plane resolution is for the event planes 
constructed from sub-events (half-events), which is most 
relevant for the present studies because the event planes 
of the half-events was used. Figure [23] shows the 6ep 
as a function of centrality. The event-plane resolution is 
maximal in mid-central collisions and decreases towards 
peripheral and central collisions. 

To study the dependence of the present results on 
the event-plane resolution, a certain fraction of particles 
were randomly discarded (/=25%, 50%, and 75%) from 
the event-plane reconstruction, thereby artificially reduc- 
ing the event-plane resolution. Figure [24] shows S(A 2 ), 
8(A+A-) and their differences between UD and LR as a 

function of e E p(/) = y/ (cos 2(V»EP,rj>o(/) - ^EP,r,<o(/))) 
for 20-30% central Au+Au collisions. The rightmost data 
point corresponds to the event-plane resolution of the 
half-events where no particles were discarded. 

The UD and LR asymmetry correlations vary with the 
EP resolution in opposite directions. The differences be- 
tween UD and LR, A(A 2 ) and A(A + A_), increase with 
the event-plane resolution as expected. The present re- 
sults are reported as measured with respect to the recon- 
structed event-plane, without correcting for the reduc- 
tions due to the event-plane resolution. The reason this 



FIG. 23: The second-harmonic event-plane resolution of sub- 
events (77 < and 77 > 0) as a function of the number of 
participants. The statistical errors are smaller than the sym- 
bol sizes. 



correction is not performed is because it is not known 
how the asymmetry correlations depend on the EP res- 
olution outside the measured range of EP resolution in 
Fig. I24I As an estimate, the asymmetry correlation dif- 
ferences between UD and LR with a perfect EP resolu- 
tion are described in Sec. IIVAI using a linear extrapola- 
tion. The linear extrapolation with fixed zero intercept 
(dashed lines in Fig. [Ml (right panel)) would be correct if 
the high-order harmonic terms in Eq. (IA2[) are negligible. 

Despite the difficulty in extrapolating to a perfect EP 
resolution, the true differences in the correlations be- 
tween UD and LR, with ideal event-plane resolution of 
unity, should be larger than the measured differences re- 
ported here. The conclusions made using the presently 
measured correlation differences between UD and LR 
can, therefore, are only made stronger if the reaction 
plane could be measured precisely. 

The particle multiplicity asymmetry correlations were 
measured as a function of the anisotropy of those par- 
ticles within one half of the TPC relative to the EP re- 
constructed from the other half of the TPC. The event- 
plane resolution varies with the particle anisotropy i>2 bs 
even though the two quantities are from different re- 
gions of phase space. The half-event used for EP re- 
construction was randomly subdivided into two quar- 
ters, a and b. The EP of the quarter events was re- 
constructed. The EP resolution of the half event was 
assessed via (cos2(V>EP,a — V"EP,b))- Figure [53] (upper- 
left panel) shows (cos2(V'EP,a — V'EP.b)) as a function of 
W2 bs - For significantly negative 7j2 bs events, the values of 
(cos2(?/;EP,a — V'EP.b)) are negative, suggesting that the 
reconstructed EP for those events is not the true reac- 
tion plane, perhaps even being orthogonal to, rather than 
aligned with, the reaction plane. This would mean that 
the UD and LR hemispheres are flipped for the events 
with significantly negative v% hs . 

The asymmetry correlations in those events with 
nearly zero v% hs were also studied. It was found that 
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FIG. 24: (Color online) Charge multiplicity asymmetry correlations, 5(A 2 ) (left panel) and 6(A+A-) (middle panel), 
and their differences between UD and LR (right panel) as a function of the event-plane resolution cep(/) = 
\J (cos 2(t/jEP,»7>o(/) — V , EP,r ! <o(/))) in 20-30% central Au+Au collisions. The error bars are statistical only. The solid lines are 
free linear fits to the data, while the dashed lines are linear fits with the intercept fixed to zero at cep(/) = 0. 



(cos2(V'EP,a — ^EP,b)} is positive for those events for all 
centralities. Figure [35] (lower left panel) shows the esti- 
mated event-plane resolution, •y/2(cos2(V'EP,a — V'EP.b)), 
for events with |t>2 bs | < 0.04 as a function of centrality. 

The asymmetry correlation dependence on the event- 
by-event particle distribution anisotropy was also stud- 
ied with respect to the first harmonic event plane re- 
constructed from the ZDC-SMD signals. The event- 
plane resolution was obtained by the correlation between 
the event planes reconstructed from the east and west 
ZDC-SMD separately. The corresponding resolutions are 
shown in Fig. [25] (right panels) . 



Appendix C: Systematic Uncertainties 

This section describes the study of the systematic un- 
certainties. The present charge asymmetries were ana- 
lyzed by rotating the reconstructed EP by 45°. It was 
found that the UD and LR asymmetry correlations are 
identical, (Afj D ) = (A£ a ) and (A+A-)vu = (A+A-) LR , 
as expected. The present charge asymmetries were also 
calculated by randomly discarding a fixed fraction of the 
particles. Essentially the same 8{A 2 ), 8{A + AJ), A(A 2 } 
and A(A + A_) results were obtained. 

To check for possible directed flow effects, the charge 
asymmetry correlations within \rj\ < 0.5 were calcu- 
lated using the event plane constructed by particles in 
0.5 < | r]\ < 1. It was found that the observed asymmetry 
correlation results were consistent with those obtained 
from — 1 < r\ < and < rj < 1. 

The greater inefficiency in two of the sectors in the east 
half of the TPC introduces larger detector effects in the 
measurement of (A±). After subtracting the statistical 
fluctuations and detector effects, the dynamical asymme- 
try variances S(A 2 t ) were consistent between the ry > 
and rj < regions, as well as between 8(A 2 ,) and 8(A 2 _) 
from each rj region (see Fig. I22|) . The average variances 
between 8(A+) and 8(A 2 _) from the two ij regions are thus 



reported, and the maximum difference of the individual 
results from the average is considered as part of the sys- 
tematic uncertainties. The covariances, S(A + A-}, were 
also consistent between the two 77 regions, so this paper 
reports the average 8(A + AJ) including half of the differ- 
ence as part of the systematic uncertainties. 

The 0- independent track reconstruction efficiency were 
not corrected for because it does not affect the present 
asymmetry measurements. Correcting for the track re- 
construction efficiency (which is a function of 77, pt, and 
centrality) does not significantly affect the values of (A 2 ), 
8(A 2 ), A(A 2 ), (A+A_), SiA+AJ), and A(A+A^). 

The present asymmetry correlations were also studied 
by varying the event and track quality cuts. Specifically, 
the event primary vertex position, Z xtv was restricted to 
within ±15 cm (default ±30 cm). The maximum dis- 
tance of closest approach, DCA, cut was also varied be- 
tween 1-3 cm (default 2 cm) and the minimum number 
of fit points, ATfit requirement was varied between 15-25 
(default 20). For these different cut sets, the correspond- 
ing ^-acceptance corrections were used. The changes in 
the present results from the different cuts are generally 
small, but are included in the systematic uncertainties. 
The data from the full magnetic field setting (FF) and 
the reversed full magnetic field setting (RFF) were also 
analyzed separately. The results are generally consistent 
within the statistical errors. The difference between the 
FF and RFF results from the combined FF and RFF data 
set, whichever is larger, is taken as part of the systematic 
uncertainties. 

Table [IT] shows the sources and magnitudes of the sys- 
tematic uncertainties on the charge multiplicity asym- 
metry correlation measurements for two selected central- 
ity bins. The systematic uncertainties from the various 
sources are added in quadrature to yield the total sys- 
tematic uncertainties. The systematic uncertainties are 
taken to be symmetric between the positive and negative 
sides. The total systematic uncertainties are shown in 
the shaded areas in Fig. [3] and Fig. [4] 

The systematic uncertainties on the charge separation 
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FIG. 25: (Color online) The TPC second-harmonic (upper left) and ZDC-SMD first harmonic (upper right) event-plane 
resolution as a function of vf B in 20-40% central Au+Au collisions. The TPC second-harmonic (lower left) and ZDC-SMD 
first harmonic (lower right) event-plane resolution for events with \v% s \ < 0.04 as a function of centrality. The error bars are 
statistical only. 



TABLE II: Sources and magnitudes of ± systematic uncertainties. All the numbers have been multiplied by the corresponding 
number of participants Apart. The upper section is for the 40-30% centrality (Apart = 78.3) and lower section is for the top 5% 
centrality (A part = 350.6). 



Source 


<5(^UD> 


5(Ain) 


<5(^1 + ^_)UD 


5{A+A-} LR 


A (A 2 ) 


A(A+A.) 


Magnetic field polarity (FF vs RFF) 


0.0004 


0.0001 


0.0041 


0.0034 


0.0005 


0.0007 


Primary vertex z? vtx cut (15 cm vs 30 cm) 


0.0006 


0.0021 


0.0048 


0.0009 


0.0015 


0.0039 


DCA cut (1 cm vs 3 cm) 


0.0006 


0.0009 


0.0013 


0.0012 


0.0015 


0.0016 


Min. number of fit points Afl t (15 vs 25) 


0.0003 


0.0017 


0.0015 


0.0001 


0.0020 


0.0014 


TPC side (West vs East) 


0.0010 


0.0028 


0.0010 


0.0009 


0.0035 


0.0001 


Total 


0.0014 


0.0040 


0.0067 


0.0038 


0.0046 


0.0045 


Magnetic field polarity (FF vs RFF) 


0.0034 


0.0040 


0.0006 


0.0015 


0.0006 


0.0021 


Primary vertex Z vtx cut (15 cm vs 30 cm) 


0.0015 


0.0050 


0.0045 


0.0005 


0.0035 


0.0040 


DCA cut (1 cm vs 3 cm) 


0.0017 


0.0032 


0.0009 


0.0007 


0.0015 


0.0016 


Min. number of fit points Nat (15 vs 25) 


0.0031 


0.0037 


0.0023 


0.0021 


0.0006 


0.0021 


TPC side (West vs East) 


0.0073 


0.0032 


0.0022 


0.0027 


0.0084 


0.0005 


Total 


0.0089 


0.0087 


0.0056 


0.0039 


0.0093 


0.0053 



parameter A and its dependence on the observed event- 
by-event second harmonic parameter t>2 bs are also stud- 
ied in the same way. The A values for all events and 
for events with |«2 bs | < 0.04, and the intercepts and 



slopes of linear fits to A(i>2 bs ) obtained using different 
cuts and etc. are compared. The differences are assigned 
as asymmetric systematic uncertainties. They are listed 
in Table U as a function of centrality. 



[1] I. Arsene et al. (BRAHMS Collaboration), [2] B. B. Back et al. (PHOBOS Collaboration), 
Nucl. Phys. A757, 1 (2005). 



20 



Nucl. Phys. A757, 28 (2005). 

J. Adams et al. (STAR Collaboration), [27] 
Nucl. Phys. A757, 102 (2005). 

K. Adcox et al. (PHENIX Collaboration), [28 
Nucl. Phys. A757, 184 (2005). [29 
T. D. Lee, Phys. Rev. D 8, 1226 (1973). 
T. D. Lee and G. C. Wick, Phys. Rev. D 9, 2291 (1974). [30 
P. D. Morley and I. A. Schmidt, Z. Phys. C26, 627 
(1985). [31 
D. Kharzeev, R. D. Pisarski, and M. H. G. Tytgat, 
Phys. Rev. Lett. 81, 512 (1998). [32 
D. Kharzeev, Phys. Lett. B 633, 260 (2006). 
D. Kharzeev and A. Zhitnitsky, Nucl. Phys. A797, 67 [33 
(2007). 

K. Fukushima, D. E. Kharzeev, and H. J. Warringa, [34 
Phys. Rev. D 78, 074033 (2008). 
D. E. Kharzeev, L. D. McLerran, and H. J. Warringa, [35 
Nucl. Phys. A803, 227 (2008). 

S. A. Voloshin, Phys. Rev. C 70, 057901 (2004). [36 
B. I. Abelev et al. (STAR Collaboration), 
Phys. Rev. Lett. 103, 251601 (2009). [37 
B. I. Abelev et al. (STAR Collaboration), 
Phys. Rev. C 81, 054908 (2010). [38 

B. Abelev et al. (ALICE Collaboration), 
Phys. Rev. Lett. 110, 012301 (2013). [39 
Quan Wang, Ph.D. thesis, Purdue University, 2012 

http: / /drupal. star. bnl.gov /STAR/theses / phd/quanwang [40 
larXiv:1205.4638| . [41 
K. H. Ackermann et al. (STAR Collaboration), Nucl. In- 
strum. Meth. A499, 624 (2003). 

F. S. Bieser et al. (STAR Collaboration), Nucl. In- [42 
strum. Meth. A499, 766 (2003). [43 

C. Adler et al., Nucl. Instrum. Meth. A499, 433 (2003). 
J. Adams et al. (STAR Collaboration), 
Phys. Rev. Lett. 92, 112301 (2004). 

B. I. Abelev et al. (STAR Collaboration), 
Phys. Rev. C 79, 034909 (2009). 

K. H. Ackermann et al. (STAR Collaboration), 
Nucl. Phys. A661, 681 (1999). 

M. Anderson et al, Nucl. Instrum. Meth. A499, 659 
(2003). 

A. M. Poskanzer and S. A. Voloshin, Phys. Rev. C 58, 
1671 (1998). 

Gang Wang, Ph.D. thesis, UCLA, 2005 



http://drupal.star.bnl.gov/STAR/theses/ph-d/gang-wang 
L. Adamczyk et al. (STAR Collaboration^ 
Phys. Rev. Lett. 108, 202301 (2012). 
F. Wang, Phys. Rev. C 81, 064902 (2010). 
S. Pratt, S. Schlichting, and S. Gavin, Phys. Rev. C 84 
024909 (2011). 

J. Adams et al. (STAR 
Phys. Rev. Lett. 95, 152301 (2005). 
M. M. Aggarwal et al. (STAR 
Phys. Rev. C 82, 024912 (2010). 
B. I. Abelev et al. (STAR 
Phys. Rev. Lett. 102, 052302 (2009). 
B. I. Abelev et al. (STAR 
Phys. Rev. C 80, 064912 (2009). 
B. I. Abelev et al. (STAR 
Phys. Rev. Lett. 105, 022301 (2010). 
H. Agakishiev et al. (STAR 



arXiv: 1010.0690 

H. Petersen, T. Renk, and S.A. Bass 
014916 (2011). 
L. Adamczyk et al 



Collaboration), 
collaboration), 
Collaboration), 
Collaboration), 
Collaboration), 
Collaboration), 
Phys. Rev. C 83, 
Collaboration), 
and B. Miiller, Phys. Rev. C 



(STAR 

larXiv: 1302.38021 
M. Asakawa, A. Majumder, 
81, 064912 (2010). 

A. Bzdak, V. Koch, and J. Liao, Phys. Rev. C 81, 
031901(R) (2010). 

G.-L. Ma and B. Zhang, Phys. Lett. B 700, 39 (2011). 
Note that in events with zero ellipticity there can be other 
harmonic shapes, such as triangularity. In other words, 
events with zero ellipticity are not necessarily spherical. 
S.A. Voloshin, Phys. Rev. Lett. 105, 172301 (2010). 
The present track reconstruction algorithm jumps over 
gaps between sectors to search for the next hits when 
reconstructing tracks across the sector boundaries. For 
track trajectories nearly parallel to a sector boundary, the 
empty distance over the gap between the adjacent sectors 
that the algorithm has to bridge is relatively large. The 
search window for hits in the next sector is made propor- 
tionally wider. In central collisions, where the TPC hit 
occupancy is high, the confusion in the hit finding greatly 
increases. This results in a lower efficiency in reconstruct- 
ing those tracks nearly parallel to the sector boundaries 
in more central collisions. 



